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The status of neutrino oscillation searches employing nuclear reactors as sources 
is reviewed. This technique, a direct continuation of the experiments that proved 
the existence of neutrinos, is today an essential tool in investigating the indications 
of oscillations found in studying neutrinos produced in the sun and in the earth's 
atmosphere. The low-energy of the reactor De makes them an ideal tool to explore 
oscillations with small mass differences and relatively large mixing angles. In the 
last several years the determination of the reactor anti-neutrino flux and spectrum 
has reached a high degree of accuracy. Hence measurements of these quantities at 
a given distance L can be readily compared with the expectation at L = 0, thus 
testing i/g disappearance. While two experiments, Chooz and Palo Verde, with 
baselines of about 1 km and thus sensitive to the neutrino mass differences associated 
with the atmospheric neutrino anomaly, have collected data and published results 
recently, an ambitious project with a baseline of more than 100 km, Kamland, is 
preparing to take data. This ultimate reactor experiment will have a sensitivity 
sufficient to explore part of the oscillation phase space relevant to solar neutrino 
scenarios. It is the only envisioned experiment with a terrestrial source of neutrinos 
capable of addressing the solar neutrino puzzle. 



a: 

Oh; 



X 



I. INTRODUCTION 

Neutrinos have the distinction of being the first el- 
ementary particle whose existence was predicted by 
a theorist in order to explain seemingly unrelated 
phenomena*. Pauli made such prediction in 1930 in his 
famous letter in order to explain the continuous electron 
energy distribution in nuclear beta decay. It became im- 
mediately clear that neutrinos will be difficult to observe, 
because the corresponding cross sections are so tiny. But 
in a series of experiments in 1953-59 Reines and Cowan 
(1953, 1959) were able to prove convincingly that electron 
anti-neutrinos from nuclear reactors are able to cause 
the inverse neutron beta decay, Pe + P —^ s^ + n, and 
hence that they are real particles. Shortly afterwards, in 
1962, the separate identity of the muon neutrinos, f^, was 
demonstrated (Danby ei a/. 1962). Another decade later, 
in 1975, the t lepton was discovered (Perl et al. 1975) 
and the observation of its decay properties implied the 



*For early development in neutrino physics see, for example, 
Chapter 1 in Winter (1981). 



existence of the third neutrino, z/,-, that was directly ob- 
served only very recently (Kodama et al. 2001). Precise 
measurements of the decay width of the Z have shown 
that just three neutrino flavors (2.994 ± 0.012 from the 
combined fit to all LEP data (Groom et al. 2000)) par- 
ticipate in the weak interactions (at least for neutrinos 
with masses less than l/2Mz). 

Phenomenologically, it is obvious that neutrinos of 
each flavor are either massless or at least many orders 
of magnitude lighter than the corresponding charged lep- 
tons with which they form the weak interaction doublets. 
Based on these empirical facts, the Standard Electroweak 
Model postulates that all neutrinos are massless, and con- 
sequently have conserved helicity (which is the same as 
chirality in this case) and that the separate lepton num- 
bers for electron, muon, and tau flavors are conserved. 
Challenging this postulate of the vanishing neutrino mass 
has recently become a central issue in many disciplines 
of fundamental science, including particle and nuclear 
physics, cosmology, and astrophysics. The present re- 
view is devoted to one particular aspect of this broad 
effort. 

Ironically, while our knowledge of intrinsic neutrino 
properties remains quite poor, these particles have been 



used as tools to understand other phenomena. The tradi- 
tion of underground neutrino detectors began thirty years 
ago when Davis and his collaborators were first able to 
detect neutrinos from the Sun. (For the description of 
the history of solar neutrino research see (Bahcall 1989).) 
This was, and together with the other experimental ob- 
servations of the solar neutrinos still is, the only clear 
proof that the basic energy generation in stars is under- 
stood. The birth of neutrino astronomy can be associ- 
ated with the observation of the neutrino burst from the 
supernova 1987A. Neutrino-induced reactions played an 
important role in establishing what is now known as the 
Standard Model of electroweak interactions when in 1973 
the neutral currents were discovered via the observation 
of the i^fj, + e -^ i^fj. + e scattering as well as the neu- 
tral current scattering of neutrinos on nucleons. Finally, 
neutrinos have been extensively used in deep inelastic 
scattering experiments at CERN and FNAL, exploring 
the quark structure of nucleons. 

The main problem in neutrino physics today is the 
question whether neutrinos, like all charged fermions, 
have a mass^. Since direct kinematic tests of neutrino 
mass lack at present the required sensitivity, the recent 
hints for neutrino mass are indirect, based on the phe- 
nomenon of neutrino oscillations. If neutrinos are mas- 
sive particles which behave in analogy to quarks, the 
states with a definite mass (i.e., the "mass eigenstates" 
which propagate as plane waves in vacuum) are not nec- 
essarily the partners of the charged leptons that couple 
to the vector bosons W^ in doublets (i.e., the weak or 
fiavor eigenstates) 
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The weak eigenstates \i^i) will be in such a case linear 
superpositions of the mass eigenstates \i>i) 
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where the coefficients Ui^i form the leptonic mixing ma- 
trix. If we assume that only three neutrinos can con- 
tribute in the Eq.(2) above, then t/ is a unitary 3x3 
matrix-'- . 

If Eq.(2) is valid, we encounter the phenomenon of neu- 
trino oscillations in which a neutrino that was initially in 



^For an up-to-date discussion of the neutrino masses and 
the relevant experiments see (Fisher, Kayser, and McFarland 
1999) and (Zuber 1998). 

^Sometimes more than 3 mass eigenstates are considered. 
The additional neutrinos must be sterile or heavy, i.e., they 
must not participate in weak interactions so that the con- 
straint from the invisible width of the Z boson is obeyed. 



the weak eigenstate I can be spontaneously transformed, 
at least in part, into another weak eigenstate neutrino of 
flavor V . (The idea of oscillations was discussed early on 
by Pontecorvo (1958, 1967) and by Maki, Nakagawa and 
Sakata (1962).) 

To see how that transformation happens, recall that 
the mass eigenstate \vi) propagates according to the ex- 
pression 

\v,{t)) ^ e-*(^-*-P'^) 1^,(0)) ~ e-^("'/2^)^|i.,(0)) , (3) 

where L is the flight path and in the last expression we 
assumed that the laboratory momenta and energies are 
much larger than the neutrino rest masses rrii. Let us 
consider now the propagation of a neutrino which was 
created at L = as a weak eigenstate \vi). At a distance 
L this state is described by 

|^.;(L))^^C/Me-("?/2^)^|^.) 
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In the last expression we used the inverse transforma- 
tion to Eq.(2), i.e., from the mass eigenstates back to 
the weak eigenstates. This step must be taken since the 
only way one can detect neutrinos is through their weak 
interactions. And in order to detect the neutrino flavor 
we have to use the charged current weak interactions, 
characterized by the production of the charged leptons 
\V). 

Thus, the neutrino of flavor / acquired components cor- 
responding to other flavors I' . This is a purely quantum 
mechanical effect, a consequence of the coherence in the 
superposition of states in Eq.(2). The probability that 
the "transition" I -^ V happens at L is obviously 
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This is an oscillating function of the distance L. The os- 
cillation length depends on the differences of the neutrino 
mass squares, \m1 — m^\, and the oscillation amplitude 
depends on the mixing matrix U . 

The necessary and sufficient conditions for the exis- 
tence of neutrino oscillations is then the nonvanishing 
value of at least one neutrino mass rrii and the nonvan- 
ishing value of at least one nondiagonal matrix element of 
the mixing matrix U. If these conditions arc hflfiUed, the 
individual lepton flavor numbers (electron, muon, tan) 
are no longer conserved. 

There is no fundamental theory which would allow us 
to deduce the parameters describing the mixing matrix U 



and the mass differences Am? . These unknown parame- 
ters must be determined empirically, by various neutrino 
oscillation experiments. Such analysis is often performed 
in a simplified way by assuming that only two neutrino 
flavors mix, e.g. e and fi. The mixing matrix U then 
depends only on one mixing angle 9, and the oscillation 
probability, Eq. (5), is also simplified 
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TTij, and we assume, as before, that 



the neutrinos are ultrarelativistic. The probability that 
fo remains Vc is obviously 



P{v,^v,,L) = l-P{v, 
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In this two-flavor scenario the oscillation amplitude is 
s\v?29 which vanishes if = or 90° and is maximum if 
9 =45°. The oscillation length is 



J^nsc — ^TT 



2E^ _ 2.48£;^(MeV) 



meters 
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To test for oscillations, one can perform either an ap- 
pearance search in which one looks for a new neutrino fla- 



vor (i.e., the deviations of P(i^o 



, L) from zero), or a 



disappearance test in which one looks for a change in the 
flux normalization (i.e., the deviation of P{i'c ^ i^c,L) 
from unity). In either case, tests performed at distance 
L are only sensitive to the values of Am^ for which 
L > 0{Losc)- Or, in other words, neutrino oscillations 
are observable only when Am'^L/E ^ 0(1). 

So far we have considered only propagation of neutri- 
nos in a vacuum. When neutrinos propagate in matter, 
such as in the solar interior, the oscillation pattern may 
be modified. This happens because electron neutrinos 
can forward scatter on electrons by charged current in- 
teractions, and other neutrino flavors cannot. Under fa- 
vorable circumstances a resonance enhancement of the 
oscillation amplitude, the so-called Mikheyev-Smirnov- 
Wolfenstein (MSW) effect (Wolfenstein 1979, 1980) and 
(Mikheyev and Smirnov 1986a, 1986b), can take place. 
Analogous matter induced oscillations can distinguish 
the hypothetical sterile neutrinos, which have no weak 
interactions at all, and the f^ or i>r neutrinos which in- 
teract with matter (electrons and quarks) by the neutral 
current weak interaction. Neither of these kinds of mat- 
ter effects is relevant for reactor neutrinos. 

For completeness, it is worthwhile to mention here two 
other issues important in the study of the neutrino intrin- 
sic properties. One of them deals with the charge con- 
jugation properties of the neutrinos. Unlike the charged 
leptons, which are Dirac particles, with distinct antiparti- 
cles, neutrinos can be either Dirac or Majorana particles. 
In the latter case of truly neutral neutrinos, there is no 



distinction between the neutrinos and their antiparticles, 
and even the total lepton number is not conserved. In 
order to decide between these two possibilities, one has 
to look for processes that violate the total lepton num- 
ber, such as the neutrinoless double beta decay. Other 
processes of this kind, like the I'e -^ i^e oscillations (e.g. 
in the present context the emission of Ve from the nu- 
clear reactor) are typically kinematically supressed, and 
their observation is unlikely in forseeable future. The dif- 
ference between the Dirac and Majorana neutrinos, while 
of fundamental importance, does not influence the results 
of the reactor oscillation searches described below. 

The other issue worth mentioning is the possibility of 
the T or CP violation in neutrino oscillations (Cabibbo 
1978 and Barger, Whisnant and Phillips 1980). In order 
to establish violation of T or CP one would have to show 
that 



P{i^e' -^ i^i) + P(yt' 
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i.e., that for example the probability of Vy, oscillating into 
Vp is different from the probability of P^ oscillating into 

Ve- 

For the usual case of three neutrino flavors, one can 
parametrize the lepton mixing matrix in terms of the 
three angles 9\ ~ ^13, 92 = 923, and 9^ = 9i2 and the CP 
violating phase 5. 
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where e.g. c\ = cos^i and s\ = sin^i, etc. 
The magnitude of the T or CP violation is character- 
ized by the differences 

Pifi ^ e) - P{y. ^ e) = -\P{ji ^ f) - P{y. ^ r)] = 

P(e -^t) - P{e ^ f) == 

-4c^siC2S2C3S3sin(5[sinAi2 +sinA23 +sinA3i] , (11) 

where A^ — [mf — m'j) x L/2E. 

Thus, the size of the effect is the same in all three 
channels, and CP violation is observable only if all three 
masses are different (i.e., nondegenerate) , and all three 
angles are nonvanishing. As will be shown below, reactor 
experiments constrain the angle 9i (or ^13), to be small 
(sin^ 2^13 < 0.1). If that mixing angle vanishes exactly, 
no CP violation is observable in the lepton sector, inde- 
pendently of the value of the CP violating phase S. To 
further improve the sensitivity to sin 29i3 iS; therefore, 
a matter of utmost importance. 



II. PHYSICS MOTIVATION OF MODERN 
EXPERIMENTS 

Numerous searches for neutrino oscillations were per- 
formed during the last two decades using nuclear reactors 
as well as particle accelerators as sources. Since most 
of them have not observed evidence for neutrino oscilla- 
tions, their results are usually based on the simplified two 
neutrino mixing scenario and presented as an "exclusion 
plot", i.e., based on them certain ranges of the parame- 
ters A?7i^ and sm^26 can be excluded from further con- 
siderations as shown in Figure 1. However, at the present 
time there are three groups of measurements that suggest 
the existence of neutrino oscillations. (And, at the same 
time, the parameter ranges suggested by them are not 
excluded.) Only these positive results will be briefly dis- 
cussed here, other experiments are listed in the Review 
of Particle Properties (Groom et al. 2000). 



A. Experimental Indications for neutrino oscillations 

The most prominent group of measurements which 
are commonly interpreted as evidence for neutrino os- 
cillations are often referred to as the "atmospheric neu- 
trino anomaly" (Kajita and Totsuka 2001). Primary cos- 
mic rays impinging on the nitrogen and oxygen nuclei 
at the top of the earth's atmosphere produce mostly 
pions, which subsequently decay via the chain tt^ -^ 
/i~P^,/i^ -^ e~T>(.v^ (and the analogous chain with 7r+ 
etc.). At sufficiently low energy, when such chains can 
fully develop, the resulting atmospheric neutrinos there- 
fore are expected to follow the v^j^ : Vg — 2 : \ ratio, 
which is essentially independent of the details of the com- 
plicated process that created them. In addition, in an 
underground detector, one can deduce the direction of 
the incoming neutrinos from the direction of the leptons 
(e and /i) created by the charged current interactions. 
Again, one is reasonably confident that this zenith angle 
distribution can be accurately predicted. If the v^ and/or 
Ve neutrinos oscillate, one expects deviations from the 
2:1 ratio mentioned above. Also, since the zenith angle 
is simply related to the neutrino path length, one expects 
deviations from the expected zenith angle dependence of 
the lepton yield. 

Both signatures of neutrino oscillations were in fact ob- 
served. The v^/ve ratio is noticeably smaller, only about 
60%, than the expected value. This result has been con- 
firmed in four detectors thus far. The anomalous zenith 
angle dependence was first observed in Kamiokande, and 
has been now confirmed, with much better statistical 
significance, by SuperKamiokande (Kajita and Totsuka 
2001, Fukuda 1998, Fukuda 2000). If these effects indeed 
signify neutrino oscillations (and we do not know of an- 
other viable explanation) then the corresponding mixing 



angle is large, s\v?29 ~ 1, and the value of the mass pa- 
rameter Am^ is in the range 10~^ - 10~^ eV^. While 
the preferred scenario at present involves v^ — > Vr oscil- 
lations, it is not clear that v^ -^ v^^ oscillations are fully 
excluded. 
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FIG. 1. Phase-space for neutrino oscillations. The existing 
limits on i/c — V/i are compared with current and future ex- 
periments and the regions obtained by interpreting the solar, 
atmospheric and LSND neutrino anomalies as due to oscilla- 
tions (some of these effects are not necessarily Uc ~ z^^ oscilla- 
tions.) The MSW mechanism is used in plotting some of the 
solar neutrino regions. The sensitivity of reactor experiments 
is the same for v^ — Vr oscillations. Limits are at 90% CL. 



The second set of measurements that can be inter- 
preted as evidence for neutrino oscillations deals with the 



"missing" solar neutrinos (Kirsten 2000, latest results in 
Fukuda 2001 and in particular in Ahmad 2001). The 
Sun produces an intense flux of electron neutrinos as a 
byproduct of the fusion reactions which generate solar 
power. It is believed that the solar structure is under- 
stood sufficiently well so that the flux and energy spec- 
trum of the neutrinos can be confidently predicted. The 
solar neutrino fluxes have been measured in six experi- 
ments so far. All of them report a deficit, i.e., the mea- 
sured flux is less than the expected one. Moreover, the 
reduction depends on the neutrino energy, inferred exper- 
imentally from the thresholds of the individual detectors. 
The only viable explanation of the deficit appears to be 
neutrino oscillation (i>e disappearance). The hypothe- 
sis that solar i/g indeed oscillate into 'active' neutrinos 
that scatter on electrons via the neutral current weak 
interaction is supported, at the 3 a level, by combining 
the pure charged current measurements of SNO (Ahmad 
2001) with the charged + neutral current measurement 
of Super-Kamiokanke (Fukuda 2001). 

By contrast to the attempts to explain the deficit by 
modification of the solar model, which are unsuccessful, 
all existing data can be simply and elegantly explained 
by invoking neutrino mass. In particular, the solution 
based on the MSW effect offers the most popular sce- 
nario. Treating the problem in the two-flavor framework 
explained above, one arrives at several isolated islands in 
the A?Ti^ - sm^29 plane. Two solutions correspond to 
Am^ w 10~^ eV^. One of them ( "small mixing angle" 
or SMA) has sin^20 wlO^^, while the other one ("large 
mixing angle" or LMA) has sin^20 > 0.5. This latter 
solution, currently giving the best fit to the data, spans 
an interval of Arn^ extending up to 10^^ eV^. The other 
possibilities have large mixing angles and Am? « 10^^ 
eV2 (LOW) or Am^ « 10" i" eV^ (vacuum). 

Finally, the only indication for oscillations involving 
man-made neutrinos comes from the LSND experiment 
which finds evidence for the i?^ — > Dc and, with more 
limited statistics, also for v^^ -^ v^, (Athanassopoulos et 
al. 1995, 1996, 1998). The former channel uses neutri- 
nos from the pion and muon decay at rest, with energies 
less than TOp/2. The latter channel uses neutrinos from 
the pion decay in flight which have somewhat higher en- 
ergies. These are appearance experiments; the observed 
signal should be absent if neutrinos do not oscillate. The 
well determined quantity is the oscillation probability, 
which has the value of about 3x10^'^. This result has 
not been independently confirmed. An analogous exper- 
iment which also uses the neutrinos from the pion and 
muon decay at rest, Karmen (Armbruster et al. 1998, 



Eitel 2000), found no evidence for the z/^ 



oscilla- 



LSND signal, they all came from measurements involv- 
ing neutrinos produced by natural sources outside of our 
control. A number of new experiments has been per- 
formed or are in various stages of planning in order to 
investigate further these tantalizing effects. Reactor ex- 
periments play an all-important role in this quest, owing 
to their unique ability to investigate very small neutrino- 
mass differences. 

Like in many other aspects of neutrino physics there 
is a fundamental difference between the past reactor os- 
cillation experiments^ and the more recent experiments 
with baselines of 1 km or more: experiments in this 
latter category are designed to further investigate, in a 
controlled environment with man-made neutrinos, par- 
ticular regions of the oscillation parameter space where 
there are indications for oscillations from other exper- 
iments. Hence the results from the new generation of 
reactor Ue detectors directly impact our understanding 
of the neutrino mixing matrix. 



B. Reactor versus accelerator-based oscillation experiments 

Nuclear reactors produce isotropically 9^ in the [3 de- 
cay of the neutron-rich fission fragments. All detectors 
optimized for oscillation searches take advantage of the 
relatively large cross-section and specific signature of the 
inverse-/3-decay reaction D^ + p —>■ n + e~^ . Such cross- 
section is shown in Figure 2 as function of the neutrino 
energy along with the neutrino flux at the reactor and the 
resulting interaction rate in a detector. We note here that 
the detection reaction has a threshold of about 1.8 MeV. 
Many of the merits and limitations of reactor-based ex- 
periments can be understood observing that the energy 
of Pe is rather low, in the few-MeV range. It directly fol- 
lows that reactor-based experiments can only be of Dg - 
disappearance type since the neutrino "beam" does not 
have sufficient energy to produce muons (or taus) and the 
neutral-current reactions of the "oscillated" P^ or P^ have 
too low cross-section and are un-distinguishable from the 
many backgrounds present. This first limitation makes 
reactor-based experiments well suited only for investi- 
gating relatively large mixing angles. In practice experi- 
ments have reported mixing sensitivities around 10% at 
large Am^ (although the proposal for a very ambitious 
experiment with sensitivity better than 2% at a particu- 
lar Ato^ will be discussed later). The second limitation 
of reactor-based oscillation searches derives from the 



tions. However, the parameter space compatible with the 
LSND signal is not fully excluded by Karmen. 

As we can see from this brief discussion, the last decade 
brought us a number of clues. With the exception of the 



^For a general discussion of sliort-baseline reactor experi- 
ments see Boehm and Vogel 1992 and Boehm 2001; for the 
individual experiments see Kwon et al. 1981, Zacek et al. 
1986, Achkar 1992, Achkar et al. 1995 and 1996, Vidyakin et 
al. 1994, Alfonin et al. 1998. 
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FIG. 2. Reactor v^. flux, inverse beta decay cross section, and 
f-e interaction spectrum at a detector based on sucii reaction. 
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FIG. 3. Neutrino Am'^ sensitivity as a function of total re- 
actor power and detector fiducial mass for detection based 
on the inverse- /3 reaction discussed in the text. The baseline 
scales with the Am^ sensitivity sought according to Eq. (8). 
The fiducial-mass X power necessary for the experiment grows 
with the square of the baseline. The past experiments are 
labelled by the name of the reactor complex used. The ap- 
proximate year of the experiment is also indicated to show 
that the increased baseline and Am^ sensivity followed more 
or less the chronological order. 



fact that the only known method of coUimating neutrino 
beams employs the Lorentz boost of the parent particles 
from which decay the neutrinos are produced. For this 
reason low energy neutrinos are generally produced over 
large solid angles, while high energy ones may come in 
relatively narrow beams. Obviously a reactor emits Ve in 
a completely isotropic way, and this, together with the 
modest interaction cross-sections available at low energy, 
makes the specific signal rates rather low. At the same 
time, however, the low energy neutrinos provide us with 
a unique opportunity to probe the lowest regions of Ato^ 
that are otherwise beyond the reach of accelerator-based 
searches. Some of these tradeoffs are well illustrated by 
Figure 3 where the Arn^ sensitivity is shown, together 
with the necessary baseline, versus the reactor power and 
detector fiducial mass for different statistical accuracies. 

Oscillation searches using reactors as sources are par- 
ticularly important today since several of the indications 
for neutrino oscillations shown in Figure 1 point to re- 
gions of the parameter space at very small Arn^ and 
nearly- full mixing. Hence two reactor-based experiments, 
Chooz and Palo Verde, were performed to investigate 
the phenomenon of atmospheric neutrinos as Ve — > Vx os- 
cillations. Such experiments, described in detail below, 
had baselines of about 1 km and fiducial masses of the or- 
der of 10 tons. For comparison, the much more complex 
accelerator-based MiNOS project between FNAL and the 
Soudan mine (Wojcicki 2001a) and analogous projects 
between CERN and Gran Sasso, Opera and Icarus 
(for a brief description, see e.g. Wojcicki 2001b), will ac- 
cess similar Am^ values with GeV-energy neutrinos and 
a baseline of the order of 1000 km. However, the 5400 
ton Minos detector and its analog at Gran Sasso will be 
able to investigate also oscillation channels not including 
v^ and reach a mixing parameter sensitivity substantially 
better than 1%. 

The reactor-based Kamland experiment, with a base- 
line larger than 100 km, will offer the unique opportunity 
of testing, with man-made neutrinos, the large-mixing- 
angle MSW solution of the solar neutrinos puzzle. In 
this case the restriction to v^ -^ i'l: -oscillations does not 
limit the interest of the experiment (since solar neutri- 
nos do certainly involve z/g ) , while its Arri^ sensitivity is 
well beyond what can be practically achieved by accel- 
erators (for comparison similar Am^ sensitivity could be 
achieved in an accelerator-based experiment with base- 
lines of order 10^ km, larger than the diameter of the 
earth) . 

Of course, the relatively lower energy of neutrinos from 
reactors pushes the optimization of reactor-based ex- 
periments to concentrate on the reduction and rejection 
of backgrounds from natural radioactivity that is, on 
the other hand, hardly an issue in accelerator-based de- 



tectors. In this respect the correlated signature of the 
inverse-/? process, the detection of the e^ and neutron, 
plays a very important role. 

While in the case of neutrinos produced by accelerators 
the experimenter has full control over the status of the 
beam, the flux of Ve. cannot be changed at will in commer- 
cial power nuclear reactors. However, in practice, typical 
reactor optimization requires a refueling shutdown every 
12 to 24 months. Such shutdowns usually last about a 
month, providing a convenient flux modulation that can 
be used to validate background subtraction methods. As 
explained in detail later, even in the case of Kamland 
that observes the neutrinos from about 70 reactor-cores, 
a substantial flux modulation is provided by the coinci- 
dence of scheduled refueling outages in the Spring and 
Fall periods, when electricity demand is lowest. 

Finally, wc remark here that the fully isotropic flux 
produced by nuclear reactors eliminates the problems re- 
lated with beam pointing that arc present in experiments 
using accelerators. While the pointing accuracy required 
in these experiments is well within the present technol- 
ogy, a fool-proof cross check of the beam-detector align- 
ment is certainly not trivial to obtain. 

In conclusion reactor-based and accelerator-based ex- 
periments offer complementary approaches to the quest 
for neutrino oscillations. It is likely that only the com- 
bined efforts on these two fronts, together with other 
studies such as the search for neutrino-less-double-/3 de- 
cay, will allow us to elucidate the problem of mixing in 
the lepton sector. 



III. REACTOR NEUTRINO SPECTRUM AND FLUX 
DETERMINATION 



Bugey (Achkar et al. 1995, 1996) reactors. Indeed such 
measurements helped gaining confidence in the reactor 
yield estimates, and, although were not the main goal of 
Chooz and Palo Verde, have been recently proposed 
(Mikaelyan 2000) for a more accurate determination of 
the mixing angle ^13 for the atmospheric neutrino region. 
Chooz could take advantage of the « 115 m distance be- 
tween the reactors for deriving weaker exclusion limits, 
that were however less affected by systematics. 
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FIG. 4. Expected positron energy spectra for no oscilla- 
tions (full line) and oscillations with parameters Am^ = 
7.2 X 10"^ eV^ and sin^ 26* = 1 at the Chooz (L ~ 1 km) 
(dashed hne) and Palo Verde (L ~ 0.8 km) (dotted hne) ex- 
periments. Adapted from Harrison, Perkins and Scott (1996). 



Since reactor-based oscillation experiments are of the 
disappearance type, the accurate determination of the 
Pe spectrum and its absolute normalization are essential 
ingredients of the measurements. We note here that for 
oscillation parameters well within the experimental sen- 
sitivity the evidence for oscillations would manifest itself 
as a deficit of events accompanied by a distortion of the 
energy spectrum as shown by the example in Figure 4. 
However, as the true value of the oscillation parameters 
moves closer to the sensitivity boundary of the experi- 
ment, the spectral shape loses power and the accuracy 
of the measurements essentially relies on the total event 
count and, hence, the knowledge of the absolute reac- 
tor flux. This last scenario also corresponds to the more 
usual case in which no oscillations are observed and an 
upper-limit is set, as well as to the case of large Am^ 
where the spectrum distortions are washed out. 

While in this section we will concern ourselves mainly 
with a-priori reactor De yield determinations, multiple- 
baseline measurements are possible and have been per- 
formed in the past at the Goesgen (Zacek et al. 1986) and 



A. Anti-neutrino production 

The determination of the De yield proceeds, schemat- 
ically, in three steps. First, the thermal power of each 
reactor core is measured, accurately and essentially con- 
tinuously. Based on such measurements, and starting 
from the initial fuel composition, the burn-up state can 
be computed as function of time. Small corrections due 
to other reactor parameters that modify the criticality of 
the core are also introduced at this time. Reactor simu- 
lation codes are often used at this stage and produce an 
accurate instantaneous fission rate for each of the rele- 
vant isotopes through the fuel cycle. In the second step 
the neutrino spectrum is derived from the fission rate. 
Finally, as the last step, the neutrino spectrum emitted 
by the reactors must be converted into an estimate of 
the experimental observable, the positron spectrum in 
the detector. Each of these steps will be explained in a 
separate subsection. 



Typical modern commercial Light Water Reactors 
(LWR) have thermal powers of the order of 3 GWth- 
This figure applies to both Pressurized Water Reactors 
(PWR) and to the less common Boiling Water Reac- 
tor (BWR) designs. In both cases the fuel is enriched 
to 2-5% in ^"^^U. Since on average each fission produces 
~ 200 MeV and ~ 6 i'e we conclude that the typical yield 
is ~ 6 X 10^" De core~^ s^^ (of course part of this flux 
will be below the detection threshold, see Figure 2). 

It is easy to understand why ~ 6 i^e are produced per 
fission. Take, as an example, the most common ^35^ 
fission, which produces two unequal fragments, and typ- 
ically two new neutrons that sustain the chain reaction, 



The existence of the 1.8 MeV threshold in the detec- 



235 



U 



Xi+X2 + 2n . 



(12) 



The mass distribution of the fragments (so-called fis- 
sion yields) is shown in Fig. 5. The lighter fragments 
have, on average, ^ ~ 94 and the heavier ones A ~ 140. 
The stable nucleus with A = 94 is 40 Zr® and the sta- 
ble A = 140 nucleus is 58Ce^^°. These two nuclei have 
together 98 protons and 136 neutrons, while the initial 
fragments, as seen from the equation above, have 92 pro- 
tons and 142 neutrons. To reach stability, therefore, on 
average 6 neutrons bound in the fragments have to P 
decay, emitting the required 6 De- 



tion process I'e + P 



automatically insures that 



only De from large Q- value, and hence short half-life, (3- 
decays are detected. Thus the observed Pe signal tracks 
closely in time the power excursions in the reactor. This 
is of some practical importance as large quantities of 
spent fuel are usually stored on-site by reactor operators. 
There is no need to track the inventory of spent fuel and 
to worry about the /3-decays of the neutron-activated re- 
actor materials which have typically a low Q-value and 
therefore long half-life products. In practice, after a few 
hours from reactor turn on/off, the detectable De flux can 
be considered saturated. 



B. Fission-rates determination 

The four isotopes whose fission is the source of virtually 
all the reactor power are ^^^U, ^^*U, ^^^Pu, and ^''^Pu. 
The fission rates deriving from their evolution during a 
typical fuel cycle in one of the Palo Verde reactors is 
shown in Figure 6 as calculated by a core simulation pro- 
gram (Miller 2000). For comparison we also show the 
evolution of ^"'"Pu and ^^^Pu that give the next to lead- 
ing contributions. The contribution of these isotopes is 
of order 0.1% or less and will not be considered further. 
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FIG. 5. Yields (in %) for ^^^U thermal neutron fission (nor- 
malized to 2 for the two fragments) 



While the total number of Pe is easy to estimate, and 
can be accurately determined given the known fission 
yields, their energy spectrum, which is of primary interest 
for the oscillation searches discussed here, requires more 
care. In particular, the commonly used neutrino detec- 
tion reaction, the inverse neutron /3 decay, has ^1.8 MeV 
threshold. Only about 1.5 Pe/fission (i.e. ^ 25%) of the 
total are above that threshold and hence can be detected, 
hence the total of ~6 Pe per fission is irrelevant. 
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FIG. 6. Time evolution of fission rates for each of the six most 
important isotopes in one of the Palo Verde reactor cores. The 
horizontal scale covers a full fuel cycle, at the end of which 
about 1/3 of the core is replaced with fresh fuel. Only the 
four most important isotopes are normally used to predict 
I'e yields. 



Each isotope produces a unique neutrino spectrum 
through the decay of its fission fragments and their 



daughters, so plutonium breeding results in a small but 
noticeable change in the emitted neutrino spectrum. 

Two types of uncertainties can be attributed to the iso- 
tope compositions described in Figure 6: errors deriving 
from uncertainties in the initial fuel composition and in 
the measurement of the plant parameters that are used 
as input to the simulation, and errors due to imperfect 
core and neutronics modeling by the simulation program 
itself. The errors intrinsic to the simulation are known to 
contribute by substantially less than 1% to the neutrino 
yield from tests in which fuel is sampled and analyzed for 
isotopic composition at the end of a fuel cycle. 
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FIG. 7. Correlation between v^, yield and the five most impor- 
tant inputs to the core simulation for a PWR. The numbers 
in the key are the slopes of the fitted lines. Note that a vari- 
ation of even 10% in any of the parameters, but power, has 
little effect on the output of the simulation. 



The correlation between the v^ yield and the plant 
parameters used as input to the simulation is shown in 
Figure 7. Apart from the obvious correlation with the 
thermal power, other parameters enter the simulation be- 
cause they affect the criticality by altering the neutron 
transport in the core (generally by the water density and 
boron absorber concentration). We see that for the pa- 
rameter with largest correlation besides power, the water 
temperature in the cold legs, an error of 10% produce an 
uncertainty of only 0.15% in the v^. yield. Of course the 
inlet temperature is known to much better that 10%. 

Economic and safety reasons provide plant operators 
with an incentive to accurately measure the thermal 
power of the reactors. Indeed, usually more than one 
method is used and the results are compared to under- 
stand the size of the uncertainties. Calorimetric methods 
(simultaneous measurement of temperature and flow-rate 



of the water outlet in the secondary cooling loop) give the 
smallest error (~ 0.6 - 0.7% at Chooz and Palo Verde) 
and are used as primary power estimate. 



C. From fission rates to the v^ spectra 

The instantaneous fission rates of the four isotopes 
235u^ 238^^ 239p^^ g^j^j 24ip^ fQm^(j above are then used 

as an input for the evaluation of the v^ spectrum. For all 
but 238y careful measurements of the /3 spectrum from 
fission by thermal neutrons were performed (Schrecken- 
bach et al. 1985, Hahn et al. 1989). These are converted 
to neutrino spectra as explained below. However, Refs. 
(Schreckenbach et al. 1985, Hahn et al. 1989) do not in- 
clude 238u 1)^0^1 undergoes only fast neutron fission and 
hence was not accessible to such measurements. 

There are, at present, several methods available to 
evaluate the De spectrum. For the De associated with 
235u^ 239p^^ g^j-^j 2^iPu, which undcrgo thermal neutron 
fission, it is customary to use a hybrid method, based 
on the conversion of the measured electron spectra as- 
sociated with the thermal neutron induced fission, into 
the Ve spectra. Clearly, the electron and Dg originate 
from the same (3 decay and share the available endpoint 
energy. For a single branch the conversion is therefore 
trivial. However, in general there are many branches, 
and many nuclei with different charges. For electron and 
antineutrino energies well above the electron mass the 
two spectra are quite similar (Schreckenbach et al. 1985) 
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(13) 



where E^ = -^/pf 



is the full electron energy, and the 



proportionality constant deviates from unity by at most 
5%. Naturally, one would like to convert the spectra more 
accurately. 

Formally, the conversion can be performed exactly as 
follows. Let n{E,Z), assumed to be a continuous func- 
tion, describe the distribution of endpoints E and nuclear 
charges Z . The electron spectrum is then 

/•OO 

Y{E,) = / AEn{E, Z)k{E, Z)p,Ee{E - E,fF{E,,Z), 

(14) 

where k(E, Z) is the spectrum normalization constant, 
Pe is the electron momentum, and F{Ei,, Z) is the Fermi 
function describing the Coulomb effect on the emitted 
electron. Provided that the electron spectrum Y{Ee) is 
measured, the endpoint distribution can be determined. 
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1 



d3 



( Y{E) 
2k{E, Z) AE^ \pEF{E, Z) 



(15) 



Once the distribution n{E, Z) is known, the Vg spec- 
trum is readily calculated by the integral analogous to 
the Eq. (14). In (Davis et al. 1979, Vogel et al. 1981) it 
was shown that such conversion procedure depends only 
very weakly on the value of Z. In practice, an empiri- 
cal relation Z{E) between the average Z and the elec- 
tron energy has been used (Schrcckenbach et al. 1985) 
(Z = 49.5 - 0.7E - OmE^ with E in MeV). 

When using the measured electron spectra the above 
expressions, involving third derivatives, are obviously im- 
practical. Instead, the integral in Eq. (14) is replaced by 
a finite sum of 30 hypothetical beta decay branches with 
branching ratios bi and equidistant endpoints Eq, 

Y{E,) = J2 hk{ElZ)5{E,,El,)p,E,{El - E,fF{E„Z), 

i 

(16) 

where S{Ee,EQ) describes the small outer radiative cor- 
rections. One can now begin with the largest value of 
-Eq (only one branch) and determine the corresponding 
branching ratio bi using the electrons of energies between 
that Eq and the next smaller one, and continue in this 
fashion step by step until the smallest Eq is reached. Pos- 
sible variations in the number and distribution of the 
endpoints Eq affects the resulting i>e spectrum not more 
than at 1% level (see Schrcckenbach et al. 1985). Having 
determined the set bi, Eq it is trivial to obtain 

YiE^J = Y. ^*^(^o: Z)EI{Eq - E,J 

i 

X [{EQ-E-,^f-mlY'^F{EQ-E^^,Z) , (17) 

where the irrelevant radiative corrections were omitted. 
This is the procedure used in deriving the neutrino spec- 
tra associated with fission of ^'^^U, ^^^Pu, and ^^^Pu 
which account for about 90% of the reactor v,. . 

The Pe spectra in Refs. (Schrcckenbach et al. 1985, 
Hahn et al. 1989, Davis et al. 1979, Vogel et al. 1981, 
Klapdor and Metzinger 1982a, 1982b) are given as tables. 
A somewhat less accurate, but easier to implement, ana- 
lytical approximation is given in Vogel and Engel (1989). 
(The fit error to the total rate is about 1.2% for ^ssy ^j^j 
only about 0.3% for ^s^Pu and 24ipu (Miller 2000).) 

For the v^ associated with the ^38^ fission one has to 
use the straightforward summation of the spectra of the 
De from all individual [3^ decays. Thus 



diV 



J2YniZ,A,t)Y,bnAE'o)P^iE,,ElZ) , (18) 



where Yn{Z, A, t) is the number of /3 decays per unit time 
of the fragment Z, A after the fissioning material has been 
exposed to neutrons for a time t and the label n char- 
acterizes each fragment. For t larger than the /3 decay 
lifetime of the fragment Z, A the quantity F„ converges 



toward the cumulative fission yield and becomes inde- 
pendent of t. Naturally, each fission fuel is characterized 
by a different set of yields Yn . 

The quantities bn^iiEg) in Eq. (18) are the branch- 
ing ratios for the ith branch with the maximal electron 
energy (endpoint energy) Eq. The branching ratios are 
normalized to unity. 

Finally, the function Pp{Ep,El, Z) is the normalized 
spectrum shape. It is usually assumed, here and in the 
conversion method explained above, that all relevant f3 
decays have the allowed shape. 
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(19) 



To use Eq. (19) is a very good approximation in practice 
and causes a totally negligible error. 

The weakness of this method is the incomplete infor- 
mation on the endpoint distribution and branching ratios 
of some fission fragments, in particular those with very 
short lifetimes and high decay energies. These 'unknown' 
decays contribute as much as 25% of the Vg at energies 
above 4 MeV. In practice, nuclear models are used to 
supplement the missing data. Examples of calculations 
based on this method are (Davis et al. 1979, Vogel et 
al. 1981, Klapdor and Metzinger 1982a, 1982b, Teng- 
blad et al. 1989). An example of an extension to lower 
De energies, where the neutron activation of the reactor 
materials plays a role, is given by Kopeikin, Mikaelyan, 
and Sinev (1997). 

While the summation method played an important role 
in the early oscillation searches, at present it is needed 
only for the description of the D^ from ^38^ fission as 
pointed out above. That component contributes only 
about 11% to the neutrino signal. We show below that 
the error associated with the summation method is less 
than 10% and hence it contributes less than 1% to the 
overall uncertainty. 

The ultimate check of the accuracy of the prediction 
outlined above consists in comparing the results in terms 
of De energy spectrum with the measurements performed 
in short baseline reactor oscillation experiments. 



D. From i/e to positrons 

Since in all reactor experiments one measures the 
positron spectra, and not directly the Ue spectra, one 
has to understand quantitatively how these are related. 
In other words, one has to know the cross section of the 
'detector' reaction 9^ +p ^ e+ + n. 

The total cross section for this reaction, neglecting 
terms of order E^/M, is given by the standard formula 
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(20) 
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where Ee^' = E,^ — {Mn — Mp) is the positron energy 
when the (smah) neutron recoil is neglected, and pi 
is the corresponding monientuni. The vector and axial- 
vector coupling constants are / = 1, g = 1.26 and 
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(21) 



where the energy independent inner radiative corrections 
are Af„„,, ~ 0.024. 

The cross section can be expressed in terms of the neu- 
tron lifetime and the phase space factor f^), ~ 1.7152 
(Wilkinson 1982) as 
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(22) 



In this way, the cross section is tied directly to the neu- 
tron lifetime, known to 0.2% (Groom et al. 2000), no 
knowledge of Gp, f / g or the Cabibbo angle 6c is in fact 
needed. 

The (small) energy-dependent outer radiative correc- 
tions to atot are given in Vogel (1984) and Fayans (1985). 
The corrections to the cross section of order E^/M , which 
arc not negligible even for the reactor energies, and the 
angular distribution of the positrons are described in Vo- 
gel and Beacom (1999). The exact threshold of the reac- 
tion is 

)2 _ M2 

-^ = 1.806 McV (23) 
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FIG. 8. Upper panel: total cross section for Pe + p — > e+ + 
n; bottom panel: {cos 9); as a function of the antineutrino 
energy. The solid line is the 0{1/Mn) result and the short- 
dashed line is the 0(1) result, Eq. 20. 



Using the results of Vogel and Beacom (1999) one can 
evaluate the total cross section as well as the quantity 
(cos^e) which characterizes the positron angular distri- 
bution essentially exactly. These quantities are shown in 
Figure 8. The high energy extension of the total and dif- 
ferential cross section has been discussed already in the 
classic paper by Llewellyn-Smith (1972). Near threshold, 
however, that treatment must be modified as shown in 
Vogel and Beacom (1999). 

The positron angular distribution, characterized by 
(cos^e) is rarely accessible. It is of interest, however, to 
consider also the angular distribution of the recoil neu- 
trons that are also detected. Since in the laboratory sys- 
tem the proton is at rest, the neutron is initially emitted 
at a forward angle restricted by 
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^2E,A ~ (A2 - m2) 
E^ 



(24) 



where A = M„ - Mp -- 1.3 MeV. The average (cos(6'„)) 
is considerably closer to unity (Vogel and Beacom 1999). 

It is often possible to localize the points where the 
positron was created and where the neutron was cap- 
tured and even though the neutron undergoes many elas- 
tic scatterings before capture, its final position maintains 
some memory of its original direction. Simulations sug- 
gest that the typical displacement of the two vertices 
is (x) ~ 1.5 cm in the organic scintillator. In fact, in 
previous reactor experiments (Zacck et al, 1986, Achkar 
1992, Zacek 1984, Achkar et al. 1995) the neutron dis- 
placement was clearly observed, in the Goesgen experi- 
ment in particular, at ~ lOcr level. The same effect was 
also observed at Palo Verde. Moreover, the single ves- 
sel Chooz experiment was able to measure the average 
neutron-positron separation and to base on it a determi- 
nation of the 9e incoming direction with an uncertainty 
of ~ 8° (ApoUonio et al. 2000). 

Given a reliable simulation of the neutron transport, 
this asymmetry allows, albeit with large errors, a di- 
rect measurement of the detector background. In case 
of 9q detection from a future Supernova, this technique 
may provide, as shown by ApoUonio et al. (2000), a crude 
but useful determination of the direction to the Vg source, 
i.e. of the star position. 



E. Accuracy of the flux and spectrum predictions 

Once the cross section and De spectra are known, the 
corresponding positron yield is easily evaluated. In reac- 
tor experiments, the neutron recoil is quite small (10-50 
keV) and thus the positron energy is simply related to 
the incoming De energy. 
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where, as before, A = M„ — Mp and we used cos6'e = 
as a good approximation of the average (cos 9e)- We note 
here that possible detector efficiency dependence on the 
positron energy requires special care. 
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FIG. 9. Positron spectrum observed by the Goesgen exper- 
iment for three different baselines (Zacek 1984). The two 
continuous lines represent fits to the data and the predictions 
obtained, as described in the text, using the measurements 
(Schreckenbach et al. 1985, Hahn et al. 1989) for ^^^U, ^^^Pu, 
and ^'*^Pu and theoretical calculations for '^''^U. 



In order to check the accuracy of the prediction one has 
to compare the results, in terms of 9^ energy spectrum 
and total flux normalization, with the measurements per- 
formed in short baseline reactor oscillation experiments. 
Since such experiments have not reported the observation 
of oscillations, we can assume that their measurements 



represent the direct determination of the reactor spec- 
trum at production. There are four factors needed for 
the evaluation of the expected positron yield and spec- 
trum: the distance to each core, the number of target 
protons, the cross section for 9^ + P -^ n 4- e+, and of 
course the quantity one wants to test, the De spectrum 
at the source. 

The distance to the reactors is trivially obtained with 
negligible uncertainty. The number of protons in the tar- 
get requires knowledge of the chemical composition and 
mass of both scintillator and possible other detector ma- 
terials where De can be captured and recorded with finite 
efficiency. Typically errors smaller than 1% are achiev- 
able for this parameter. The cross section has been dis- 
cussed at the previous subsection, and its uncertainty is 
also less than 1%. Finally, the De spectrum is the object 
of the test. 

The De flux normalization alone was tested in Dcclais 
et al. (1994) where the total neutron yield was measured 
with high statistical accuracy, and found to be in agree- 
ment with the expectation based on the known neutron 
lifetime at the 3% level. 

The validity of the tests performed at short distances 
from a reactor is reinforced by the fact that some of 
the experiments such as Goesgen (Zacek et al. 1986) 
or Bugey (Achkar et al. 1995, 1996) did measurements 
at different baselines observing no difference between 
the spectra. This is shown for Goesgen in Figure 9. 
Moreover, the relatively recent Bugey 3 measurements 
were performed at 15-40 m distance from the core and 
recorded very high statistics (some 1.5 x 10^ De events). 
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FIG. 10. Ratio between Bugey 
ent predictions. In a) the measurements are compared to 
the a-priori calculations of Klapdor and Metzinger (1982a, 
1982b). In b) Bugey 3 data is compared to the prediction 
obtained using the (3 spectra measurements of Schreckenbach 
et al. (1985) and Hahn et al. (1989), and the calculation 
mentioned for ^'^^U. The dashed envelopes are estimates of 
the overall systematics. Adapted from Achkar et al. (1996). 
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The good agreement between Bugey 3 data and the 
non-oscination predictions is demonstrated in Figure 10. 
In panel a) the prediction is generated purely from the- 
ory. Of more practical importance is panel b) where the 
prediction derives from (3 spectra (except for ^ssy -^^i^gj-e 
theory is used). In this case a fit to a horizontal line gives 
a level of 0.99 with ^jd.o.j. = 9.2/11. The final Bugey 3 
result is quoted as having a 1.4% total uncertainty. 
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FIG. 11. Expected number of Vf, interactions in the Palo 
Verde detector during the ~2 years of data-taking of the 
experiment. Note that one reactor is closer to the experiment 
while the other two are equidistant; this explains the different 
excursion for one of the refuelings. The steady decline in 
Ve interactions during the cycle is the effect of fuel burn-up. 



For reference we give in Figure 11 the time evolution 
of the r?e interaction rate expected in the Palo Verde 
detector, calculated as described above from the plant 
data. This time evolution is typical of a plant with more 
than one reactor. Refueling outages (about 1 month long 
each) give different rate excursions due to the different 
distances between the reactors and the detector. Short 
accidental reactor trips are also visible along with the 
steady rate decline through the reactor cycle due to fuel 
burn-up. 

In Fig. 12 the measured effect of the changing reactor 
fuel composition is shown. At the same time, the figure 
demonstrates how small that effect really is. We plot in 
Fig. 12 af/Ef since the number of events Ui, at a given 
time and fuel composition is 



1 Wth 
47ri?2 {Ef) 



Np£(Tf , 



(26) 



where R is the distance, Np the number of protons, e 
is efficiency for the event detection, c/ is the effective 



cross section per fission and {Ef) is the average energy 
per fission which are both sensitive to the burn-up. 

In conclusion, the v^ spectra, and its absolute normal- 
ization, are known to about 2% accuracy. Obviously, the 
reactors as Ve sources are perfectly isotropic. The differ- 
ences between various reactors, and the time changes due 
to the fuel burn-up are small, and well understood. They 
do not cause additional uncertainty. Thus, for the reactor 
neutrino oscillation searches at a few percent accuracy, 
no short distance 'monitor' detectors are needed. One 
can simply compare the measured positron spectra at 
the distance L to the expectation at the source {L — Q). 
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FIG. 12. The measured v^ event rate at Chooz, expressed as 
the cross section per fission divided by the energy per fission 
(see text for explanation), as a function of the reactor burn- up 
(in units of MWxdays per ton), i.e., the accumulated reac- 
tor power per unit mass of fuel. The fits to the expected 
rate, which includes the reactor burn-up, and to the fiat rate, 
where the burn-up is neglected, are also shown. Adapted from 
Nicolo (1999). 



F. Other reactor neutrino experiments 

We shall only briefly mention here other experiments 
with reactor v^ that are of importance as a potential 
source of information about neutrino intrinsic properties. 
Besides protons, two targets have been extensively stud- 
ied: deuterons and electrons. 

Reactor v^ can cause deutcron disintegration by two 
reaction channels; the charged current (CC) 

Pe + rf^e++n + n, (27) 

with the threshold of 4.03 MeV (M„ + me - Mp + B), 
where B is the deutcron binding energy, and the neutral 
current (NC) 
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p 



(28) 



with lower threshold oi B = 2.23 MeV. The cross sec- 
tions of these reactions, ^ 1.1 x 10^'*^ cm^/fission for the 
CC, and ^ 3.1 x 10^^^ cm^/fission for the NC (see e.g. 
Davis et al. (1979) and Vogel et al. (1981)) are more 
than an order of magnitude smaller than for the reaction 
on the proton target. (The cross section is expressed 
'per fission' because the fission rate is the quantity more 
directly related to the reactor power than the v^ fiux.) 

The study of the Ve + d reactions was pioneered by 
Reines and collaborators (Pasierb et al. 1979, Reines et 
al. 1980), who observed the corresponding two or one 
neutron captures. From the point of view of neutrino os- 
cillations, the ratio of the CC and NC rates is potentially 
useful, since the CC is flavor sensitive and NC is not (the 
same idea is being pursued in the case of solar neutri- 
nos by the SNO collaboration) . More recent experiments 
(Riley et al. 1999, Kozlov et al. 2000) show, as expected 
given the short distance from the reactor, no indication 
of oscillations. 

The other reaction observed with the reactor v^ is the 
scattering on electrons 



e , 



(29) 



where the spectrum of the recoil electrons (originally as- 
sumed at rest in the laboratory since the electron mo- 
mentum associated with the atomic binding is usually 
negligible) is observed. Obviously, the reaction signature, 
just the recoiling electron, is quite difficult to distinguish 
from background caused by radioactivity, making the ob- 
servation of the Ue — e scattering very challenging. 

The cross section for v^ — e scattering consists of the 
well understood weak interaction part, and a so far un- 
observed incoherent electromagnetic part: 



dcr 
dT 
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(30) 



where T is the kinetic energy of the recoiling electron. 
The first part, weak scattering, represents the sum of co- 
herent (interfering) contributions from charged and neu- 
tral currents, while the second part, proportional to a^, 
can be described as due to a finite neutrino magnetic 
moment /i,y. Only massive neutrinos can have magnetic 
moments, and hence the study of this reaction, and the 
possible determination of /i^, is of great interest. 

Again, pioneering results were obtained by Reines and 
collaborators (Reines, Gurr and Sobel 1976). More re- 
cently, limits on /^^ with reactor neutrinos, albeit with a 
rather poor signal to noise ratio, were obtained by Derbin 
et al. (1993). At present a slightly more stringent direct 



limit of fi^, < 1.6 X 10~^°y^B comes not from the reac- 
tor neutrinos, but from the analysis of the shape of the 
Super-Kamiokande (Fukuda et al. 1999) solar neutrino 
data (Beacom and Vogel 1999). A new effort to improve 
the sensitivity to /z^ is currently underway (Amsler et al. 
1997) at the Bugey reactor (MUNU experiment). 



IV. EXPERIMENTS MOTIVATED BY THE 
ATMOSPHERIC NEUTRINO ANOMALY 

Two experiments have been built with the specific pur- 
pose of testing the hypothesis that neutrino oscillations 
occur with the parameters found by the atmospheric 
neutrino measurements. Both experiment are now com- 
pleted. All Chooz data are published in ApoUonio et 
al. 1998, 1999, and Nicolo 1999. Data-taking at Palo 
Verde finished in the summer 2000 and the results are 
published in Boehm et al. 2000a, 2000b, 2001, and Miller 
2000. 

As it can be seen from Figure 3, in order to access 
Am? « 10~^ eV^ with reactor neutrinos, a baseline of 
order 1 km and a mass in excess of a few tons are needed. 
Indeed, the backgrounds from cosmic radiation and natu- 
ral radioactivity are a major consideration in the design 
of such large, low-energy detectors, and different back- 
grounds situations led the two groups to rather different 
designs. 

The Chooz detector was built in a pre-existing un- 
derground cavity under a ~ 100 m rock overburden 
(~ 300 m.w.e). This substantial cosmic radiation shield- 
ing allowed the use of a homogeneous detector where 
inverse-/? events were tagged as a double (delayed) co- 
incidence between the e"*" (prompt) signal and the n (de- 
layed) one. Such simple event signature can be identified 
with large efficiency so that a 5 ton active mass was suffi- 
cient for the experiment. The Palo Verde detector, on 
the other hand, was located in an underground bunker 
excavated for the purpose. Economic considerations lim- 
ited the overburden to 12 m (~ 32 m.w.e.) sufficient only 
to eliminate the hadronic component of the cosmic radi- 
ation and reduce the muon fiux by a factor of five. The 
rather large remaining muon fiux produced a substan- 
tial quantity of secondary neutrons so that a segmented 
detector was needed to take full advantage of the triple 
coincidence given by the e"*" ionization and subsequent 
7's from annihilation. This more elaborate topological 
signature reduced the detector efficiency and pushed the 
fiducial mass to 12 ton. Both detectors were built with 
materials selected for low radioactivity and included a 
passive 7 and neutron shield and an active cosmic-ray 
veto counter. 

While the 3-reactor plant of Palo Verde produces a 
larger power (11.6 GWth) than the 2-reactor one at 
Chooz (8.5 GWt^), more important is the fact that, 
unlike in the case of Palo Verde which was a plant 
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already running for a long time, the Chooz reactors 
were commissioned after the start of data taking of the 
experiment. This endowed the collaboration with the 
rather unique opportunity of observing the backgrounds 
at reactor-off for a substantial period and the slow ramp- 
up of power. On the other hand the need to cope with a 
much more stable operation, with the periodic ^ 2/3 
refueling power excursions shown in Figure 11, moti- 
vated the Palo Verde group to develop new methods 
for background subtraction that will be important for 
future experiments likely to run in similar steady-state 
situations. 
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FIG. 13. Initial light attenuation length in Palo Verde scin- 
tillator. Top: attenuation length as function of wavelength in 
three pre-production batches. Bottom: attenuation length at 
A = 440 nm for the initial batches of scintillator. Produc- 
tion batches are to the right of the vertical dashed line, while 
on the left we show various test batches not used in the de- 
tector. The first point ("CF") refers to a standard (non Gd 
loaded) scintillator fiuid. Typical production scintillator had 
an attenuation length better than 11m. 



is not easy to achieve in the stable and sufficiently trans- 
parent form needed for a large detector. Both groups in- 
vested substantial resources in scintillator development. 
In Figure 13 we show initial attenuation length data for 
the Palo Verde scintillator that was a cocktail of 60% 
mineral oil, 36% pseudocumene (1,2,4 trimethylbenzene) 
and 4% alcohol (used to keep the Gd compound in solu- 
tion), with PPO as primary fluor. This scintillator had a 
H:C ratio of ~ 2 and a light yield of 56% of anthracene, 
with typical light attenuation length greater than 11 m 
at A = 440 nm (Piepke et al. 1999). 

The time stability of the same scintillator is shown 
in Figure 14 where the light attenuation curve for one 
cell and the effective attenuation length for all cells are 
presented as measured three times, at one year intervals. 
The 12% loss in the first year decreases to 3% in the 
second, possibly indicating a gradual stabilization. 
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FIG. 14. Evolution of scintillation-light attenuation-length in 
Palo Verde scintillator during the two years of life of the 
detector. Top: the attenuation curve for scintillation events 
at different locations along one of the 9 m-long cells. Bottom: 
effective attenuation values for all detector cells. Note that 
the shorter value of the effective attenuation length reflects 
the non-trivial optics of the cells. 



Both detectors used liquid scintillator loaded with 
0.1% natural gadolinium which has a high thermal neu- 
tron capture cross-section and releases a large amount 
of energy in the capture. In this way the neutron cap- 
ture time is reduced to ^27 fj,s from ~170 fj,s for the 
unloaded scintillator, proportionally reducing the uncor- 
related background. Furthermore, Gd de-excitation after 
the capture releases a 8 MeV 7 cascade, whose summed 
energy gives a robust event tag well above natural ra- 
dioactivity. In contrast, neutron capture on protons re- 
sults only in a single 2.2 MeV 7. While the Gd loading 
offers obvious advantages in suppressing backgrounds, it 



The time variation of the attenuation length Acd of the 
Gd-loaded scintillator at Chooz, regularly measured in 
the detector, is well fitted by the empirical function 



AGd(t) 



Ao 



1 + at 



(31) 



which accounts for the observed exponential decrease of 
the signal with time. Results of the fit are shown in 
Figure 15. 
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The properties of the Gd-loaded and unloaded veto 
liquid scintillators developed for Chooz are presented in 
Tabic I. 





A,(3/t) = yd + at) 
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FIG. 15. Light attenuation length Xcd of the Chooz Gd- 
loaded scintillator versus time with best-fit function shown. 



TABLE I. Main properties of the liquid scintillators used in 
the Chooz experiment. 





Gd-loaded 


unloaded 


Chem. cont.: 
basic 

aromatics, 
alcohols 
wavelength 
shifters 


Norpar-15** 

(50% vol.) 

IPB-l-hexanol 

(50% vol.) 

p-PTP+bis-MSB 

(1 g/1) 


Mineral oil 
(92.8% vol.) 

IPB 

(7.2% vol.) 

PPO + DPA 

(1.5 g/1) 


Atomic mass 
composition: 
H 
C 
Gd 
others 


12.2% 
84.4% 
0.1% 
3.3% 


13.3% 

85.5% 

1.2% 


compatibility 


acrylic, tefion 


density (20°) 


0.846 g/1 


0.854 g/1 


Flash point 


69° 


110° 


Scint. yield 


5300 photons/MeV (~ 35% of anthracene) 


Optical 
atten. length 


4 m 


10 m 


Refr. index 


1.472 


1.476 


Neutron 
capture time 


30.5/is 


180^s 


Thermal neutron 
capture path 


~ 6 cm 


~ 40 cm 


Capture 
fraction on Gd 


84.1% 





A. Chooz 

The Chooz detector was built at the distances of 1115 
m and 998 m from the two reactors of the new Chooz 
power plant of Electricite de France in the Ardenne region 
of France. The plant, shown in Figure 20, has a total 
thermal power of 8.5 GWt^ and the two reactors reached 
full power in, respectively. May and August 1997. The 
experiment took data from April 1997 until July 1998, in 
the conditions specified in Table II. 

The apparatus, schematically shown in Figure 16, con- 
sisted of a central volume of scintillator with a mass of 
5 tons, where Dg were detected. This scintillator was 
contained in an acrylic vessel (region 1) that separated 
it from a 70 cm thick shielding layer of mineral oil (re- 
gion 2). 192 eight-inch photomultiplicrs (PMT's) were 
mounted onto a steel vessel that, in-turn, isolated, me- 
chanically and optically, the central detector from the 
outer veto counter. The central detector had a photo- 
cathode coverage of 15% and a light yield of ^ 130 pho- 
toelectrons/MeV (p.e./MeV) (Baldini et al. 1996). The 
90 ton veto scintillator was at least 80 cm thick and 
was readout with two rings of 24 eight-inch PMT's, the 
outer containment tank being painted with white reflec- 
tive paint. An outer layer (75 cm thick) of low- activity 
sand provided primary shielding from the rock. 

Laser flashers were installed to monitor the detector 
performance and radioactive sources could be inserted 
into the central region of the detector through special 
pipes. The detector energy response was calibrated daily 
with ^°Co, ^^^Cf and AmBe 7- and n-sources in order to 
accurately track the aging of the scintillator, the detector 
efficiency and the energy calibration. 

As an example we show in Figure 17 the results of a Cf 
calibration run with the source placed in the middle of 
the detector. Data is compared with a Monte Carlo sim- 
ulation for the reconstruction of the x, y and z positions 
and total energy in the detector. Both the peaks for n- 
captures on p (2.2 MeV) and Gd (8 MeV) are clearly vis- 
ible. The very good energy resolution {a{E)/E — 5.6 % 
at 8 McV) allows one to verify that the 8 McV peak is 
in fact the superposition of a 7.77 MeV line with 77% 
weighting from capture on ^'"'^'Gd (energy shifted from 

TABLE II. Summary of the Chooz data-taking conditions. 







Time (h) 


r 


Wthdt {GWhtu)\ 


Total run 




8761.7 








Live 




8209.3 








Dead 




552.4 








Reactor 1 ON 


only 


2058.0 






8295 


Reactor 2 ON 


only 


1187.8 






4136 


Both reactors ON 


1543.1 






8841 


Both reactors OFF 


3420.4 










*Norpar-15 is the trademark of Exxon Mobil Corporation. 
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7.94 MeV because of scintillator saturation effects) and 
a 8.31 MeV line with 23% weighting from capture on 
155 Gd (energy shifted from 8.54 MeV). The fit to these 
two Gaussians gives x^/d.o.f. = 67.6/55 while the fit to 
a single Gaussian is very poor with x^/d.o.f. = 875/58. 
The position resolution was found to be Ux — <Jy ^ <Jz ^ 
17.5 cm. 
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FIG. 16. Schematic drawing of the Chooz detector. 
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FIG. 17. Visible energy and position reconstructed for a cal- 
ibration ^^^Cf source placed inside the Chooz detector. 



B. Palo Verde 

The Palo Verde experiment was built at the Palo 
Verde Nuclear Generating Station, the largest nuclear 
plant in the Americas, ~ 80 km west of Phoenix, in the 
Arizona desert. The total thermal power from three iden- 
tical pressurized water reactors is 11.6 GWth- Two of the 
reactors were located 890 m from the detector, while the 
third was at 750 m. The shallow underground bunker 
housing the detector is shown in Figure 21 at the time of 
construction. In total 350.5 days of v^ data were collected 
at Palo Verde in the period between October 1998 and 
July 2000, covering four scheduled refueling outages as 
indicated in Figure 11. Of these, 242.2 days were at full 
power, while 21.8 days (86.5 days) had the reactor at 
750 m (890 m) off. Such data was complemented by fre- 
quent calibration runs. 

The fiducial mass, segmented for active background re- 
jection, consisted of 66 acrylic tanks filled with 0.1% Gd- 
loaded scintillator and arranged as shown in Figure 18. 
Each cell was 9 m long, with a 12.7 x 25.4 cm^ cross 
section, and it was viewed by two 5-inch photomultiplier 
tubes, one at each end. A v,, is identified by space- and 
time-correlated e"*" and n signals. Positrons deposit their 
energy in a scintillator cell and annihilate, yielding two 
511 keV 7's that, in general, will be detected in differ- 
ent cells, giving a triple coincidence. Neutrons thermal- 
ize and are captured in Gd, giving a 7-ray shower of 
8 MeV total energy also detected in more than one cell. 
The central detector was surrounded by a 1 m-thick wa- 
ter shield to moderate background neutrons produced by 
muons outside the detector and to absorb 7's from the 
laboratory walls. Outside of the water tanks were 32 
large liquid scintillator counters and two end-caps to veto 
cosmic muons. The rate of cosmic muons was approxi- 
mately 2 kHz. The pattern of muons traveling through 
veto counters and their timing relative to the central de- 
tector hits were recorded for subsequent off-line analy- 
sis. The central detector was equipped with a system of 
tubes that allowed the insertion of calibration sources in 
the small spaces between cells. In addition, a set of blue 
LEDs and optical fibers was used to produce flashes of 
light inside each of the cells. In order to reduce natural 
radioactivity, all building materials for the detector were 
carefully selected, including the aggregate (marble) used 
in the concrete of the underground laboratory. 

Both the positron and the neutron were triggered by 
a triple coincidence requiring at least one cell above a 
"high" threshold set at ~600 keV (positron ionization or 
neutron capture cascade), and two cells above a "low" 
threshold set at ~40 keV (Compton scattering from an- 
nihilation photons or neutron capture cascade tails). The 
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triple coincidences were required to be within 3x5 ma- 
trices of cells anywhere in the detector as recognized by 
a custom-made trigger processor (Gratta et al. 1997). 
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FIG. 18. Schematic view of the Palo Verde neutrino detec- 
tor. 



The efficiency calibration was based upon a primary 
measurement performed a few times per year with a cal- 
ibrated ^^Na e"^ source and an Am-Be neutron source. 
The ^^Na source mimicked the effects of the positron from 
the Ug interaction by providing annihilation radiation and 
a 1.275 MeV photon which simulated the e^ ionization in 
the scintillator. The source was placed at 18 positions in 
the detector deemed to be representative of different con- 
ditions. The results of this procedure were then rescaled 
to the e^ case using a Monte Carlo simulation. The neu- 
tron detection efficiency was measured by scanning the 
detector with the Am-Be source where the 4.4 MeV 7 
associated with the neutron emission was tagged with a 
miniaturized Nal(Tl) counter. Other calibrations, used 
to measure the detector energy response, were performed 
using the Compton edges from ^'^^Cs, ^^Zn and ^^^Th 
sources. The same Th source was also used more fre- 
quently to track the scintillator transparency, as already 
shown in Figure 14. Weekly runs of the fiber-optic and 
LED flasher systems were used, respectively, to monitor 
the gain and linearity of photomultipliers and the tim- 
ing/position relationship along the cells. 

Since the energy deposition of the 511 keV 7's in one 
cell has a sharp falling spectrum (Compton scattering) 
it was vital to have the lowest possible "low" thresholds 
in the trigger and to understand the behavior of such 



thresholds with great accuracy. This second task was 
complicated by the fact that the trigger used voltage am- 
plitudes, while only charge from integrating ADCs was 
available off-line. For this reason the detector simula- 
tion included a detailed description of the signal devel- 
opment in time. This code correctly described the shape 
of pulses taking into account scintillator light yield, at- 
tenuation length and de-excitation time; photomultiplier 
rise- and fall-time and gain; and event position along a 
cell. The simulation of the detector response to the ^^Na 
source is shown in Figure 19 and correctly describes the 
40 keV (600 keV) threshold position to within 1.4 keV 
(2.6 keV), resulting in uncertainties on the positron and 
neutron efficiencies of 4% and, respectively, 3%. 
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FIG. 19. A comparison of the trigger thresholds at Palo 
Verde from data and Monte Carlo. The data were taken 
with a ^^Na source at the center of each cell. The top por- 
tion shows the efficiency of the trigger thresholds (low and 
high) for a typical cell as a function of energy deposited; the 
bottom shows the energy at 50% efficiency for low and high 
thresholds in all 66 cells. 



C. Backgrounds 

There are generally two types of background affect- 
ing long baseline reactor experiments where the signal 
is based on the correlated e~^,n signature: uncorre- 
lated hits from cosmic-rays and natural radioactivity and 
correlated ones from cosmic-/i-induced neutrons. The 
first type can be measured by studying the time dif- 
ference between positron-like and neutron-like parts of 
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an event. More insidious are cosmic- /i- induced neutrons 
that present the same time and space correlation between 
prompt and delayed parts of the event as in Ve ■ Such 
events are schematically shown in the Palo Verde de- 
tector, in Figure 22. Neutrons are produced by cosmic-/i 
spallation and capture on the materials outside the veto 
counter. Both production mechanisms can result in ei- 
ther neutron thermalization and capture, where the ther- 
malization process fakes the prompt triple coincidence, 
or secondary neutrons production, where one of the cap- 
tures fakes the prompt triple coincidence. Conceptually 
the same situation holds for Chooz (and Kamland, as it 
will be discussed later), although differences in overbur- 
den and the simpler scheme of coincidence, numerically 
change the relative importance of different backgrounds. 
It is useful to point out that direct neutrons from the 
reactors have a totally negligible effect at the distances 
discussed here. 

Both experiments pre-select i>e candidates by requir- 
ing an appropriate topology (in space and time) for the 
prompt and delayed parts of each event and their rel- 
ative position. Such cuts insure that the spatial and 
temporal extents of the events are compatible with the 
i>e hypothesis and that events are well contained and 
measured in the detector. A general classification in 
terms of signal and different backgrounds can be conve- 
niently done by studying the correlation between prompt 
and delayed energy in Chooz for such pre-selected sam- 
ple, as shown in Figure 23. The region marked "B" in 
the Figure contains cosmic-ray muons stopping in the 
detector after entering undetected by the veto counter. 
Both prompt energy (muon ionization) and delayed en- 
ergy (Michel electron) are large. Indeed events in region 
"B" show a fast time correlation between prompt and 
delayed part, consistent with the muon lifetime. Region 
"C" is populated by the muon-spallation events discussed 
above: large prompt energy deposit from proton recoils 
in neutron thermalization is accompanied by a fix 8 MeV 
energy deposit characteristic of neutron capture. Regions 
"A" and "D" are populated by random coincidences of 
natural radioactivity hits, sometimes including a high- 
energy proton recoil from neutron scattering in the de- 
layed part (region "A"). Neutrino candidates populate 
the region framed by the darker line, as it can be seen by 
comparison between the scatter plots with reactors ON 
and OFF. 

The time elapsed between the prompt and delayed 
parts of the events is shown in Figure 24 for the Palo 
Verde data. We note that the process of n-capturc in 
the segmented detector requires the sum of two exponen- 
tials to fit properly. This is due to the fact that a fraction 
of the neutrons stop, after thermalization, in passive ma- 
terials (mainly acrylic for Palo Verde) where there is 
no Gd and the capture is a slower process. While the 
Monte Carlo gives a good fit with two exponentials, for 



data a third exponential, with longer time constant, is 
needed in the fit. Such exponential accounts for events 
initiated by uncorrelated background, having the delayed 
part triggered by cosmic rays, crossing the detector with 
a 2 kHz rate. Timing cuts are applied by both experi- 
ments to insure that events are consistent with a neutron 
capture. In addition events are rejected for a period of 
time following tracks detected in the veto counters. This 
last cut is particularly important at Palo Verde where 
the cosmic-ray rate is high. 

The availability for Chooz of data at zero power and 
with the reactors ramping up provides an independent 
way to check the magnitude of signal and background. 
The fitting procedure proceeds as follows. For each run 
the predicted number of neutrino candidates resulting 
from the sum of a signal term, linearly depending on 
the reactor effective power W* , and the background, 
assumed to be constant and independent of power ^t; so. 



N,^{B + W*X)M„ 



(32) 



where the index i labels the run number, At^ is the cor- 
responding live time, B is the background rate and X is 
the positron yield per unit power averaged over the two 
reactors. 

The results are listed in Table III, for three data taking 
periods corresponding to threshold readjustments. The 
data are also shown in a compact form in Figure 25. 

A simple subtraction of the e"*" spectra with reactor 
ON and OFF gives for Chooz the spectrum shown in 
Figure 26. The comparison of the observed distribution 
with the expected one for no-oscillations already shows 
very good agreement. 

TABLE III. Summary of the likelihood fit parameters for the 
three data taking periods at Chooz . 



period 


1 


2 


3 


starting date 


97/4/7 


97/7/30 


98/1/12 


runs 


579 -^ 1074 


1082 -> 1775 


1778 -> 2567 


live time (h) 


1831.3 


2938.8 


3268.4 


reactor-off time (h) 


38.9 


539.5 


2737.2 


WAt (GWh) 


7798 


10636 


2838 


B (counts/d) 


1.25 ±0.6 


1.22 ±0.21 


2.2 ±0.14 


X (counts/d GW) 


2.60 ±0.17 


2.60 ±0.09 


2.51 ±0.17 


x'/dof 


136/117 


135/154 


168/184 


N„ (counts/d) 
(Ofull power) 


24.8 ±1.6 


24.8 ±0.9 


24.0 ±1.6 



^^The "effective" power W is a fictitious thermal power 
corresponding to both reactors located at the reactor 1 site, 
and thus providing 9.55 GW at full operating conditions and 
at starting of reactor operation. 
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FIG. 20. Aerial view of the Chooz power plant. The detector is located in a tunnel under the hills on the bottom right of the 
photograph. 




FIG. 21. The Palo Verde underground laboratory at the time of construction (Fall 1996). 
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FIG. 22. Schematic view of two types of cosmic-/^-induced baclcgrounds and a signal event (far left) in the Palo Verde 
detector. Neutrons are produced by cosmic-/^ spallation (right) and capture (left) on the materials outside the veto counter. 
Both can result in either neutron thermalization and capture (left), where the thermalization process fakes the prompt triple 
coincidence, or secondary neutrons production (right), where one of the captures fakes the prompt triple coincidence. 
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FIG. 23. Delayed energy vs. prompt energy in pre-selected Chooz events. The selection cuts are listed in the figure. On the 
left it is shown the case of reactors OFF, while on the right is the case of reactors ON. A description of the event-types in the 
different regions of the plot is given in the text. 



TABLE IV. Summary of results from the Palo Verde experiment (Boehm et al. 200f). Uncertainties are statistical only. 
N, N' and (1 — £i)-Bpn are measured rates, while Background and R^ are efficiency corrected (assuming that the background 
events are measured with the same efficiency as the signal). -Rcaic is for no oscillation hypothesis. See text for notation. 



Period 


1998 


1999-1 


1999-11 


2000 


Reactor 


on 


890 m off 


on 


750 m off 


on 


890 m off 


on 


890 m off 


Duration (d) 


30.4 


29.4 


68.2 


21.8 


60.4 


29.6 


83.2 


27.5 


efficiency (%) 


8.0 


8.0 


11.5 


11.6 


11.6 


11.6 


10.9 


10.8 


N id-') 


39.6 ± f .f 


34.8 ± I.f 


54.9 ± 0.9 


45.1 ± 1.4 


54.2 ± 0.9 


49.4 ± 1.3 


52.9 ± 0.8 


43.1 ± 1.3 


N' (d-i) 


25.f ± 0.9 


21.8 ± 0.9 


33.4 ± 0.7 


32.0 ± 1.2 


32.5 ± 0.7 


32.6 ± 1.0 


30.2 ± 0.6 


30.4 ± 1.1 


(I - ei)Bp„ (d-i) 


0.88 


0.89 


1.11 


1.11 


1.11 


1.11 


1.07 


1.07 


Background (d"'^) 


292 ± f f 


255 ± fO 


265 ± 6 


266 ± 10 


256 ±6 


265 ± 9 


249 ± 5 


272 ± 9 


R. (d-i) 


202 ± 19 


182 ± 18 


212 ± 10 


124 ± 17 


214 ± 11 


161 ± 15 


237 ± 10 


129 ± 16 


i?calc (d-^) 


216 


154 


218 


129 


220 


155 


218 


154 
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FIG. 24. Time elapsed between the prompt and delayed parts 
of events in Palo Verde data and Monte Carlo. The sim- 
ulated data are fit to two exponentials. Real data are fit to 
three exponentials of which the third accounts for the random 
background. 
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FIG. 25. Chooz Ue rate during the reactors commissioning. 
The background at reactors off is f .1 ± 0.25 events/day. 



FIG. 26. Background-subtracted positron energy spectrum in 
Chooz data. Error bars represent statistical errors only. The 
solid histogram represents the expectation for the case of no 
oscillations. The ratio between the two curves is shown on 
the bottom panel. 



The same procedure can be repeated for Palo Verde 
using the thermal power excursions due to refueling. 
However in this case this technique substantially mag- 
nifies the errors since: 1) the periods of low power still 
have about 2/3 of the full flux, so that in the subtraction 
most of the signal is lost, 2) the statistical errors are dom- 
inated by the relatively short periods of low power. In 
addition, for any experiment, the background subtraction 
method will give correct result only if special attention is 
paid to the data quality, guaranteeing in particular that 
the efficiencies for signal and background are accurately 
known and remain as constant as possible through the 
experiment. 

An alternative method (Wang et al. 2000) was de- 
veloped for the Palo Verde analysis starting from the 
evidence that, for their depth and detector configura- 
tion, the dominant correlated background has at least 
two neutrons, each triggering the detector with its cap- 
ture. Such intrinsic symmetry can be used to cancel 
most of the background directly from data and com- 
pute the remaining components from Monte Carlo sim- 
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ulations. This technique makes the best possible use 
of the statistical power of all data collected. The rate 
of candidate events after all cuts can be written as 
N — i?unc + Bnn + Spn + Si, where the contribution of 
the uncorrelated i?unc, two-neutron B^n and other corre- 
lated backgrounds Bp^ are explicitly represented, along 
with the Pe signal S^. The dominant background B^^ 
(along with -Bunc) is symmetric under exchange of sub- 
events, so that an event selection with the requirements 
for the prompt and delayed event parts swapped, will re- 
sult in a rate N' = Bunc + ^nn + ei^pn + i2S^. Here 
ei and 62 account for the different efficiency for selecting 
asymmetric events after the swap. One can then calcu- 
late A^—A^' = (1 — ei)i?pn-l-(l — e2)'S'iy where the efficiency 
correction £2 — 0.2 can be estimated from the i>c Monte 
Carlo simulation. 

The Palo Verde group found that the processes of ^- 
spallation in the laboratory walls and capture of the ^'s 
that are not tagged by the veto counter, contribute to 
(1 — ei)_Bpn, while other backgrounds are negligible. Us- 
ing Monte Carlo simulation, they obtain (1 — ei)i?pn = 
-0.9 ± 0.5 d-i (-1.3 ± 0.6 d"^) for ^-spallation in the 
1998 (1999) data-set; the same figures for /i-capture are 
0.6 ±0.3 d-i (0.9 ±0.5 d^^) in 1998 (1999). This repre- 
sents only a small correction to N ~ N' since the error on 
Bpn is reduced by the fact that ei is close to 1. While the 
Monte Carlo model is accurate for the capture process, in 
the case of spallation the broad range of spectral indexes 
for the n-rccoil energy reported in literature was simu- 
lated (Wang et al. 2000). The average between different 
predictions is then used for Bp^ while the spread is used 
as an extra systematic error. Since no Vc signal is present 
above 10 MeV, the observed integrated rate above such 
energy is used as a normalization of the Monte Carlo. 
The rate of neutrons produced by muon spallation has 
been measured at Palo Verde (Boehm et al. 2000c) 
and the dependence of the neutron spallation yield on 
depth has been analyzed by Wang et al. (2001). 

The Palo Verde results obtained in this way are 
shown in Table IV for different running periods. Clearly 
also in this case there is good agreement with the no- 
oscillation hypothesis. 



best use of the PMT charge and time information. The 
relative importance of the time and charge information 
in optimizing the detector spatial resolution (which also 
affects the energy resolution) depends on several factors: 
the size of the vessel, the distance between PMTs and 
the fiducial target volume, the scintillator light yield and 
attenuation length. In Chooz, the relatively small de- 
tector volume and the small distance between the neu- 
trino target and the PMTs, made the charge information 
the dominant one in the precise event characterization. 
In the case of larger volume detectors, filled with non 
Gd loaded scintillator, like Kamland (Alivisatos et al. 
1998), or also, like the solar neutrino experiment BOREX- 
INO (Ahmonti et al. 1998, 2000), the time information 
gains importance. 

As an example of a minimization algorithm for event 
reconstruction, based on the PMT measured charge, we 
describe the one used for Chooz. The standard algo- 
rithm uses a maximum likelihood method to reconstruct 
the energy E and the vertex x of an event. The like- 
lihood is defined as the joint Poissonian probability of 
observing a measured distribution of photoelectrons over 
24 "patches", each grouping 8 adjacent PMTs, for given 
{E, x) coordinates in the detector. So, for an event oc- 
curring at time t after the start of data taking, one can 
build a likelihood function as follows: 

24 24 -JrNj _ 

C{N-N) = Y{P{N,-N,{E,x,t)) = \{^L^e "^ 



i=i 



i=i 



N,\ 



(33) 



where Nj is the observed number of photoelectrons and 
N j the expected one for the j-th patch given an event 
{E,x,t). The reason for using a Poissonian instead of 
Gaussian statistics is due to the frequent occurrence of 
low energy events with low number of photoelectrons de- 
tected by some PMT patches. The predicted number of 
photoelectrons for each patch is computed by considering 
a local deposit of energy, resulting in a number of visi- 
ble photons which are tracked to each PMT through the 
different attenuating Region 1 (Gd-doped) and Region 2 
scintillators. Therefore 



D. Event reconstruction techniques 

The identification of the neutrino signal and the rejec- 
tion of the background in reactor neutrino experiments 
depends on the accuracy of the event energy and posi- 
tion determinations and on the spatial and time corre- 
lations of the detected positron and neutron. Event re- 
construction in a segmented detector like Palo Verde 
is relatively simple. However, event characterization in 
single vessel detectors, like Chooz or Kamland, which 
are viewed by photomultipliers placed at the vessel sur- 
face, requires refined reconstruction methods making the 



N, 



aErjY^ 



fe=i 



An 



■ exp 



d^jx) dj'^ix) 



(34) 



where E is the ionization energy deposited in the scin- 
tillators, a is the light yield of the scintillator, 77 is the 
average PMT quantum cfhciency, ^j^ is the solid angle 
subtended by the k-th PMT from the event position, dj 
is the light path length in region 1, cij is the light path 
length in region 2, Xcd is the attenuation length in region 

1 scintillator, and Xui is the attenuation length in region 

2 scintillator. 

The solid angle is approximated by the expression 
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rijfe ^ 111 I 1 



Ijk 



(35) 



jk 



^PMT ^^^ 



rpMT being the PMT photocathode radius, djk — d-[ + 
d2 , and 6 being the angle between the event-PMT di- 
rection and the inward unit vector normal to the PMT 
surface. 

Instead of (33), as is usually the case for problems in- 
volving the maximum likelihood method, it is more con- 
venient to use the theorem on the "likelihood ratio test" 
for goodness-of-fit to convert the likelihood function into 
the form of a general x^ statistics (Eadie et al. 1971). 
If one assumes Nj to be the best estimate of the true 
(unknown) photoelectron distribution and can form the 
likelihood ratio A defined by 



A: 



C{N-N) 
C{N;N) 



(36) 



The "likelihood ratio test" theorem states that the "Pois- 
sonian" x^, defined by 
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x' = -21ogA = 2^[iV, 



N,+N,logi^)l (37) 



asymptotically obeys a chi-square distribution (Baker 
and Cousins 1984). It is easy to prove that the minimiza- 
tion of x^ is equivalent to maximization of the likelihood 
function, so that the x^ statistic may be useful both for 
estimating the event characteristics and for goodness-of- 
fit testing. 

The CERN-MINUIT package (James 1994) was used 
to minimize (37). The starting value for the i-th coordi- 
nate was based on the charge asymmetries measured by 
initially grouping the PMTs into only 6 "superpatches" , 
referred to the detector frame axes; it was defined ac- 
cording to: 



given for Chooz should probably be considered as some 
sort of ultimate limit for reactor-based oscillation experi- 
ments, at least when only one detector is present. Indeed 
the intrinsically high efficiency (~ 70%) of the homoge- 
neous detector, together with the unique opportunity of 
studying the zero power case, are important advantages 
(for comparison the efficiency of the larger but segmented 
Palo Verde detector is ~ 11%). 

The (energy averaged) ratio between j/g detected and 
expected was found to be 

R = 1.01 ± 2.8%(stat) ± 2.7%(syst) Chooz (39) 

and 
R = 1.01 ± 2.4%(stat) ± 5.3%(syst) Palo Verde 

(40) 

in both cases consistent with unity. 

TABLE V. Origin and magnitude of systematic errors in 
Palo Verde and Chooz. Note that the two experiments 
offer different breakdowns of their systematics. For simplicity 
we do not show the systematics for the Palo Verde ON-OFF 
analysis. The Palo Verde results are from the analysis of 
the full data set (Boehm et al. 2001). 



Systematic 


Chooz (%) 


Palo Verde (%) 


cr{i^c + P — > nH-e+) 


1.9 




- 




Number of p in target 


0.8 




- 




Wth 


0.7 




- 




Energy absorbed per fission 


0.6 




- 




Total rate prediction 




2.3 




2.1 


e+ trigger eff. 


- 




2.0 




n trigger cfi'. 


- 




2.1 




Ue selection cuts 


- 




2.1 




(1 — ei)i3pn estimate 


- 




3.3 




Total Ue efficiency 




1.5 




4.9 


Total 




2.7 




5.3 



XiO 




1,2,3, 



(38) 



where the indices H — refer to the opposite superpatches 
of the i-th axis and D' is the half size of the detector along 
that axis. Once the Xio corresponding to the starting 
position is known, the starting energy value is obtained 
from (34) after replacing x with xq and Nj with Nj. 
Examples of the results obtained by this procedure are 
shown in Fig. 17. 



E. Results and systematics 

A summary of systematic errors for both Chooz and 
Palo Verde is given in Table V. The systematic error 
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Both experiments were able to exclude v^ - v^ os- 
cillations as dominating for the atmospheric neutrino 
anomaly. This is evident from the exclusion contours ob- 
tained using the unified approach (Feldman and Cousins 
1998) and shown in Figure 27 for Chooz and Figure 28 
for Palo Verde. 



two flavor scenario, assuming that rn^ ^ 




analysis A 



— analysis B 



- - analysis C 

Q 90% CL Kamiokande (mulli-GeV) 
H 90% CL Kamiokande (sub+multi-GeV) 



0.6 0.7 



sin"(2e) 



FIG. 27. Limits on mass difference and mixing angfe from 
Chooz (90% CL) obtained witli tlie unified approach (Feld- 
man and Cousins 1998). Analysis A refers to the curve 
obtained by a fit to the background-subtracted spectrum 
in which both shape and normalization are used. Analy- 
sis C uses only the shape of the spectrum. Finally analy- 
sis B uses the difference of baselines between the two reactors 
(AL = 116.7 m). While in this last case most systematics can- 
cel, statistical errors are larger and the Am^ sensitivity rather 
poor due to the short baseline difference. The Kamiokande 
fc — Vfi atmospheric neutrino result is also shown. 



F. Are smaller mixing angles within experimental reach ? 



Amfg 



Am^g 



(i.e. 



Am , while ATO22 — 0). This sce- 
nario is obviously compatible with the evidence based 
on the atmospheric neutrino anomaly (A?n^ '^ 3 x 10"'^ 
cV^) and the solar neutrino deficit {Am? < 10"** eV^). 
In such a case the mixing angle 612 becomes irrele- 
vant and one is left with only three unknown quanti- 
ties: A?Ti^, ^13, and ^23- With this parameterization, and 
assuming that p behave like D, the Ve disappearance is 
governed by 



Pii^e -^ Vx) = sin^ 26*13 sin' 



4E, 



(41) 



while v^ — *■ Vr oscillations, in this scenario responsible 
for the atmospheric neutrino results, are described by 



P(t'p ^ i^r) = cos"* 6*13 sin^ 26*23 sin^ 



Ato^L 
4.E, 



(42) 
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FIG. 28. Limits on mass difference and mixing angle from 
Palo Verde at 90% CL (Boehm et al. 2001). The full curve 
is obtained with the swap background subtraction method de- 
scribed in the text, while the dashed is obtained using the re- 
actor power changes to estimate and subtract the background. 
The Kamiokande v^ — Vf. atmospheric neutrino result is also 
shown for illustration. 



The current data on neutrino oscillations suggests the 
need to include at least three neutrino flavors when 
studying results from experiments. As discussed in the 
Introduction, the most general approach would involve 
five unknown parameters, three mixing angles and two 
independent mass differences. However, an intermedi- 
ate approach consists of a simple generalization of the 



An analysis of the atmospheric neutrino data based on 
these assumptions has been performed (Okamura 1999) 
and its results arc shown in Figure 29 for the Ve disap- 
pearance channel. One can see that, while the relevant 
region of the mass difference Am^ is determined by the 
atmospheric neutrino data, the mixing angle ^13 is not 
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constrained very much. Here the reactor-based neutrino 
oscillation experiments play a decisive role. 

The determination of the angle 613 has obvious im- 
portance not only for the structure of the lepton mixing 
matrix U but for the observability of CP violation in 
the lepton sector, as stressed in the Introduction. If ^13 
vanishes, or is very small, no CP violation effects are ob- 
servable in the lepton sector. Moreover, for the vanishing 
^13 and with three neutrinos only, the lepton mixing is 
radically simplified. The electron neutrino is then simply 



cosfyi2i^i +Sm6'i2l'2 



(43) 



while the i/^ and i>. 



neutrinos become superpositions of 
1^3 and the corresponding orthogonal combination of vi 
and 1^2 ■ It is therefore interesting to ask whether reactor- 
based experiments can be extended to address regions of 
even smaller mixing parameter sin 2^13. 
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FIG. 29. Exclusion plot showing the allowed region of 613 and 
Am? based on the Super-Kamiokande preliminary analysis 
(the region inside the dotted curve). The region excluded 
by the neutrino reactor experiments are to the right of the 
corresponding thick and thin continuous curves. 



A simple inspection of Table V shows that, using the 
Chooz systematics, if all flux and cross-sections related 
errors could be set to zero one would be left with an error 
of ~ 1.5%. Hence, assuming a detector large enough to 
produce negligible statistical error the total error would 
shrink from the present 3.9% to 1.5%. 

This scenario is considered by Mikaelyan (2000) that 
proposes to use an underground reactor at Krasnoyarsk 
in Russia as a source and two identical detectors placed 
at distances of ~ 1100 m and ~ 250 m. The interesting 
feature of the Krasnoyarsk site is that there are substan- 
tial facilities available underground, with an overburden 
of ~ 600 m.w.e., twice the depth of Chooz. Indeed 



it might even be conceivable to locate the detectors on 
rail-cars and periodically switch their position to further 
reduce some of the systematics related to detector ef- 
ficiency. The proposal discusses the use of 50 tons of 
Gd-loaded scintillator for each of the two identical ho- 
mogeneous detectors, so that the far detector would col- 
lect 50 events/day (the thermal power of the reactor is 
in this case lower than at Chooz or Palo Verde). The 
background is estimated by Mikaelyan (2000) to be of 
5 events/day or less. 

This proposal estimates that such an experiment could 
reach a sensitivity in mixing strength of better than 0.02 
in the Am^ region relevant for atmospheric neutrinos. 
While the idea looks certainly interesting, it would be 
useful to explore how practical it is in general to push 
the errors of the absolute De flux to the 1% domain, 
even with the measurements considered here. Note also 
that the Krasnoyarsk reactors, according to the Gore- 
ChernomyrdinW agreement, are supposed to be shutdown 
for re-coring in not very distant future. 



V. EXPLORING THE SOLAR 1^ ANOMALY ON EARTH: 
KAMLAND 

While historically solar neutrinos provided the first 
hint for oscillations, there is a consensus today that the 
strongest evidence for oscillation is the atmospheric neu- 
trino anomaly. Indeed, the zenith angle dependence of 
the anomaly has substantially helped to eliminate expla- 
nations not based on some property of neutrinos them- 
selves and the advent of K2K, to be followed soon by 
the Minos and CERN to Gran Sasso programs (Wojci- 
cki 2001a, 2001b), are bringing the study of oscillations 
in this regime to a laboratory activity with both source 
and detector well under control. 

In the case of solar neutrinos none of the effects that 
would be generally considered "smoking guns" for oscil- 
lations has yet clearly emerged from the data and their 
exploration "in a laboratory setting" is made particularly 
challenging by the huge L/E^, required. It is probably a 
safe prediction that it will take a very long time before 
an accelerator-based experiment will be able to tackle the 
solar neutrino problem! However, the very low energy of 
reactor neutrinos make a reactor-based oscillation experi- 
ment able to reach the Large Mixing Angle (LMA MSW) 
solution possible - albeit rather challenging. While the 
analysis of current and future solar neutrino experiments 
presumably will help to decide which of the 



"** "US-Russian Plutonium Production Agreement is Signed" 
Statement by the White House office of the Vice President, 
23 Sept. 1997. 
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FIG. 30. Location of nuclear power plants in the World. Substantial concentrations of reactors are in Europe, East of the US 
and Japan. (Note that the map, see http://www.insc.anl.gov/, contains few plants that were either planned but never built or 
are no longer operational.) 
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FIG. 31. Location of large nuclear power plants in Japan, Korea and Far East Russia. (See the comment in the caption of Fig. 
30.) 
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possible solutions is the right one, we find the chance to 
study solar neutrino oscillations "in the lab" extremely 
compelling. Note that, unlike in the case of atmospheric 
neutrinos, where it turned-out that electron neutrinos 
are not involved in the dominant mixing, the solar neu- 
trino problem, if due to oscillations, obviously involves 
Vf. disappearance. So, unless D^ behave drastically dif- 
ferently from Vf. (which would be a worthwhile discovery 
anyway, signaling the breakdown of CPT symmetry) a 
reactor experiment is an exact replica of the astrophysi- 
cal experiment, only built on earth. 



A. Nuclear reactors in Japan 

The "easier" solution of the solar neutrino problem 
(MSW LMA) is shown in Figure 1. In order to com- 
pletely explore such solution one needs a Ato^ sensitiv- 
ity of at least 10~^ eV^ at large mixing angle. As now 
customary we refer to Figure 3 as a first step in design- 
ing our experiment: we see that a « 100 km baseline is 
needed and this drives the power x fiducial-mass product 
in the 10* MWjhXtons range. Clearly a large detector 
has to be used in conjunction with very many nuclear re- 
actors. A cursory look at the placement of nuclear power 
plants on the earth. Figure 30, reveals that such an ex- 
periment could only be placed in Europe, East of the 
United States, or Japan. 

There are 16 commercial nuclear power plants in 
Japan, their location being shown in Figure 31. They 
supply about 1/3 (or 130 GWth) of the total electric 
power in the country. At the Kamioka site there is 
an anti-neutrino flux of ~ 4 x lO^cm^^s^^ (or ~ 1.3 x 
10*^cm-2s-i for E^^ > l.SMeV) from these reactors. 80% 
of such flux derives from reactors at a distance between 
140 km and 210 km, so that there is a limited range of 
baselines. The breakdown of this data by power plant 
(several plants have on site more than one reactor) is 
given in Figure 32. We note that some 2% of the flux de- 
rives from power plants in South Korea (the Primorskaya 
plant in Russia is only planned) that will have to be in- 
cluded (albeit only as a crude estimate) to provide an 
exact flux prediction. 

While the Figure 32 assumes the nominal power for 
each of the cores, an average over one year, taking into 
account scheduled and unscheduled down times, gives an 
expected non-oscillation rate of ~750 kton~^year~^ for 
a C„H2„+2 target. Although the signal is provided by a 
very large number of cores it turns out that a modula- 
tion of the De flux is expected at Kamland (Alivisatos et 
al. 1998) thanks to the refueling and maintenance sched- 
ule of nuclear power plants in Japan. Such shutdowns, 
in fact, are concentrated in the Fall and Spring when 
the power demand is lowest, as illustrated in Figure 33. 
Hence, from the point of view of the tools available to 



study backgrounds, Kamland is in a situation very sim- 
ilar to that of Palo Verde, with 2 dips in the flux from 
full to « 2/3 expected every year. 



Site 


Distance 
(km) 


#of 
cores 


r^uier.; 
(GW) 


flux 
(v cm-^s"'') 


Signal 
(v/yr) 


Japan 
Kashiwazaki 


160.0 


7 


24, 


4.25x10^ 


348.1 


Ohi 


179.5 


4 


13.7 


1.88X10^ 


154.0 


Takahama 


190.6 


4 


10.2 


1.24x10^ 


101.8 


Hamaoka 


214.0 


4 


10.6 


1.03x105 


84.1 


Tsuruga 


138.6 


2 


4.5 


1.03x10= 


84.7 


Shiga 


80.6 


1 


1.6 


1.08x10^ 


88.8 


Mihama 


145.4 


3 


4.9 


1.03x10= 


84.5 


Fukushima-1 


344.0 


6 


14.2 


5.3x10" 


43.5 


Fukushima-2 


344.0 


4 


13.2 
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FIG. 32. List of relevant parameters for power reactors in 
Japan and South Korea. 
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FIG. 33. Power-flux level at Kamioka from Japanese reac- 
tors as function of time. Low power periods in the Fall and 
Spring seasons are alternated with peaks of high power in the 
Summer and Winter. 



It is interesting to remark that other artificial sources 
of low energy D^ are not a serious background in Kam- 
land. The largest effect would be produced by a large 
nuclear powered vessel stationed, while running its reac- 
tors at full power, in the Toyama bay, 50 - 100 km from 
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the detector. In these circumstances the excess signal in 
Kamland would amount to ~10% (Detwiler 2000). It 
is clear that it is extremely unlikely that such conditions 
will occur for any significant period of time. 



B. Detector design 

The Kamland detector is housed in the cavity built 
for the Kamiokande detector under the summit of Mt. 
Ikcnoyama in the Japanese Alps, about 50 km east of the 
town of Toyama. The layout of the laboratory is shown in 
Figure 34. The rock overburden is more than 1,000 m in 
any direction with an average rock density of 2.7 g/cm"^. 
The site is at 500 m distance from SuperKamiokande. 

A cutout view of the Kamland detector is shown in 
Figure 35. The fiducial volume consists of a sphere con- 
taining 1000 tons of liquid scintillator. The scintillator 
container is a thin plastic- walled balloon of 6.5 m radius 
that is not supposed to take the weight of the scintilla- 
tor but only to isolate it from an outer 2.5 m thick layer 
of non-scintillating, radiation shielding, fluid. The bal- 
loon is also designed to be impermeable to radon that 
mainly originates from Th and U contaminations inside 
the PMT's glass. The buffer fluid and the liquid scintilla- 
tor are contained and mechanically supported by a stain- 
less steel spherical vessel that also provides the mechan- 
ical structure where the photomultipliers for the fiducial 
volume are mounted. The sphere is solidly anchored in- 
side the cylindrical rock cavity and the space between 
them is filled with water and used as a veto Cerenkov 
counter. The scintillator, based on mineral oil and pseu- 
documene, is designed to achieve sufficient light yield and 
n-7 discrimination by pulse-shape analysis, yet comply- 
ing with rather strict flammability requirements from the 
Kamioka mine. Given the cost and stability issues for a 
detector of the size of Kamland, it was chosen not to 
Gd-load the scintillator. As it will be discussed later, 
simulations indicate that sufficient signal-to-noise ratio 
will be achieved with unloaded scintillator. Events will 
be localized inside the fiducial volume using the light in- 
tensity and propagation delays to the different photomul- 
tipliers so that large area, fast tubes are required. While 
the veto counter will be read-out using 20-inch photo- 
multipliers dismounted from the Kamiokande detector, 
new, faster, tubes with 17- inch active photocathode have 
been developed for Kamland in order to allow for proper 
vertex reconstruction from timing. Such tubes have an 
average transit-time-spread of ~ 3 ns (to be compared 
to ~ 5 ns for the Kamiokande/SuperKamiokande tubes). 
The central detector has a 30% photocathode coverage 
obtained using about 1280 seventeen-inch tubes comple- 
mented, for energy measurements, by 642 twenty-inch 
Kamiokande tubes. A spherical shell of acrylic panels 
(not shown in Figure 35) is mounted at a radius imme- 
diately inside the position of the PMT's and is used as 



primary barrier against radon migration into the active 
scintillator. A cylindrical stainless steel chimney of 3 m 
diameter protrudes from the top of the sphere to permit 
access to the central detector during installation. Buffer 
fluid and scintillator lines as well as calibration ports are 
mounted in the chimney along with all the electrical ca- 
bling. 

The readout of Kamland is designed to provide wave- 
form analysis information for each of the PMT's in the 
detector with essentially no dead-time for several consec- 
utive events. This allows for clean event reconstruction 
and enables the off-line study of the pre-history of in- 
teresting events. For example multiple neutron events, 
described above as the most dangerous background at 
Palo Verde, will be fully reconstructed by Kamland. 
Similarly cosmogenic activation giving short half-life nu- 
clei will be clearly recorded. Deep digital buffering will 
allow the detector to sustain substantial burst of events 
like expected from supernovae. 

In Figure 36 we show a phase of the central detector 
PMT installation that was concluded in September 2000. 
Scintillator filling started in Spring 2001 and data taking 
at Kamland is scheduled to begin before the end of 2001. 



C. Expected performance 

Similarly to previous experiments, both random hits 
from natural radioactivity and correlated events from 
neutron production in cosmic-ray-muon spallation and 
capture, contribute to the background to reactor Pg in 
Kamland. The results of Monte Carlo full detector sim- 
ulation using the measured Kamioka cosmic ray flux and 
the activities of various components as sampled during 
construction are given in Table VI. For the purpose of 
this background estimate U and Th contaminations in 
the scintillator of 10~^^ g/g has been assumed. Such 
purity level has been already achieved in samples of the 
Kamland scintillator. Monte Carlo studies have shown 
that cosmogenic activation gives negligible contribution 
to the background for doubles. A discussion of back- 
grounds to single signatures, not considered here, can be 
found in Alivisatos et al. (1998). 

TABLE VI. Summary of background rates in Kamland for 
the i/e signature. A signal-to-noise ratio of about 10/1 is ex- 
pected for reactor i/e ■ Adapted from Alivisatos et al. (1998). 



Background source 


Rate (day^^) 


Cosmic muons induces neutrons 

Natural radioactivity (random coincidence) 

Natural radioactivity (correlated) 


0.1 
0.15 
0.005 


Total predicted background 
Reactor Ue signal (no oscillation) 


0.25 
2 
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FIG. 34. Partial view of the system of tunnels inside Mount Ikenoyama with the locations of Kamland and its main services. 
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FIG. 35. Schematic cross-section of the Kamland detector. 



30 




FIG. 36. View of the internal volume of the Kamland sphere during the central detector installation. A modular styrofoam 
raft is used as a platform for workers. The installation begun from the top of the sphere and moved down, as the water level in 
the sphere was reduced. PMTs, black shades, acrylic plates, monitoring LEDs and cables were mounted in place for each level 
before lowering the water. 
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In Figure 37 we show the predicted energy spectra for 
reactor neutrinos at Kamland for no oscillations and 
different oscillation parameters consistent with the LMA 
MSW solar neutrino solution. We can use one of such 
curves and add to it fluctuations consistent with a 10/1 
signal-to-noise ratio and three years of data to investigate 
the sensitivity of the experiment. Assuming that oscilla- 
tions with Ato^ = 2 X 10"^ eV^ and sin^ 29 = 0.75 are 
indeed the cause of the solar neutrino anomaly, we obtain 
the measurement of the oscillation parameters shown in 
Figure 38. On the other hand, no evidence for oscillation 
after three years of data would result in the exclusion 
curve shown in Figure 1 and would rule-out the LMA 
MSW solution to the solar neutrino problem. 
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FIG. 37. Positron energy spectra expected at Kamland for 
no oscillations and oscillations with indicated parameters 
Am^ and sin^ 26 = 0.75 in the MSW LMA solar neutrino 
solution. 



D. Other physics with a very large low-energy Ue detector 

Kamland will be the largest detector specifically op- 
timized to detect low-energy v^ with good efficiency and 
low background. This opens a number of interesting op- 
portunities beyond the measurement of oscillations from 
reactors. In addition, such a large detector with a low 
energy threshold can be used to directly measure neu- 
trinos from the sun, assuming that backgrounds can be 
sufficiently reduced and understood to enable the detec- 
tion of single energy deposits. Of particular importance 
is the ^Be line that is below the threshold of the water 
Cerenkov detectors. The presence of large amounts of 
carbon in Kamland 's scintillator opens the possibility 



of detailed flavor studies in neutrinos coming from su- 
pernovae. Finally Kamland represents such a large step 
in size and backgrounds relative to the previous detectors 
that one should be ready for the possibility that it will 
discover completely new and un-expected phenomena in 
physics or astrophysics. 
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FIG. 38. Simulated measurement of neutrino oscillation pa- 
rameters from three years of data at Kamland assuming 
Am? = 2 X lO"'^ eV^ and sin^ 26 = 0.75. A signal-to-noise 
ratio of 10/1 was assumed (see text). 



Here we will only mention the topic of terrestrial 
Ve that is somewhat unusual and directly relates to the 
experiment's ability to detect v,, . The reader interested 
in the direct detection of solar neutrinos or neutrinos 
from supernovae is referred to the Kamland design re- 
port (Alivisatos et al. 1998). A description of the new 
and un-expected phenomena mentioned above will be 
hopefully provided at a later stage. 

Although the study of terrestrial anti-neutrinos was 
proposed as early as 1966 (Eders 1966) practical difficul- 
ties, due to the very small cross-sections and very low en- 
ergies involved, have made this physics impractical until 
now. Kamland has the ability to detect energy deposi- 
tions of the order of 1 MeV in a unprecedented amount of 
liquid scintillator and is therefore ideally suited for this 
study. It is important to realize that low energy v^ are 
easily detected with very low background in Kamland 
thanks to their very specific signature. 

The cooling rate of our planet and its contents of heavy 
elements are central issues in the earth sciences. The 
earth radiates about 40 TW of heat from its surface. 
About 40% of this energy (or 16 TW) is believed to have 
radiogenic origin with 90% of it deriving from decays of 
^^^U and ^'^^Th. Radiogenic heat is therefore an essential 
component of the present Earth dynamics. As discussed 
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by several authors (Eders 1966, Marx 1969, Marx and 
Lux 1970, Avilez et al. 1981, Krauss 1984) the concen- 
tration of these isotopes can be mapped, at planetary 
scale, by direct detection of Pe deriving from the /3-decay 
processes. Since neutrinos have a mean free path many 
orders of magnitude larger that the size of the Earth, 
the neutrino field is analogous to a gravitational field, 
where the sources are represented by radioactive density 
(as opposed to mass density). 

Since the maximum energy carried by terrestrial neu- 
trinos is (Krauss 1984) 3.27 MeV and the capture thresh- 
old is 1.8 MeV, the maximum in the energy spectrum de- 
tected in the prompt part of the events will be 2.49 MeV 
(including the 1.02 MeV from positron annihilations). 
For energies above threshold only the Thorium and Ura- 
nium decay chains give a detectable amount of events. 
^^''Pa from the U chain and ^^*Ac and ^^^Bi of the 
Th chain have similar endpoints (respectively 2.29 MeV, 
2.08 MeV and 2.25 MeV) while ^"Bi from the U chain 
has an endpoint of 3.27 MeV. Therefore the energy spec- 
trum observed for the prompt part of the event should 
have a characteristic double-hump structure shown in 
Figure 39. This will also allow the measurement of the 
U/Th ratio. Anti-neutrinos from nuclear reactors give, as 
described above, a similar signature, but their energy is 
substantially higher and, as shown in the Figure, they can 
be easily separated from the terrestrial anti-neutrinos. 
Indeed, a repetition of the analysis for the reactor neu- 
trinos discussed above using only positron energies above 
2.7 MeV gives an oscillation sensitivity very similar to the 
one presented in Figure 38. 
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FIG. 39. Energy spectrum from terrestrial anti-neutrinos 
compared with reactor signal as expected in Kamland. Three 
different geophysical models are shown for the terrestrial anti- 
neutrinos and no oscillations are assumed for all the spectra 
shown. Adapted from Raghavan et al. (1998). 



The two lower spectra (la and lb) superimposed in Fig- 



ure 39 for the terrestrial anti-neutrino component corre- 
spond to two different possible geophysical models with 
different heavy elements concentration in the oceanic and 
continental crusts (Raghavan et al. 1998). The highest 
curve (Ila) is given as a reference and shows what the 
spectrum would be in the extreme case where the entire 
40 TW of heat escaping from the Earth's interior was 
generated by the Th and U decay chains. 

In one year of data-taking, model la would give an 
integral of 61 events while model lb would give only 41 
events, and a differentiation between the two at 3a level 
could be obtained in five years of data-taking, taking 
into account the fluctuations of the background due to 
the reactor neutrinos. 



VI. CONCLUSIONS 

The use of nuclear reactors to study neutrino prop- 
erties has a long and glorious history. While the first 
experiments devoted to oscillation searches were moti- 
vated by the generic principle to "look where the light 
is" , many of the modern hints for neutrino oscillations 
point to parameters that match very well the capabil- 
ities of reactor-based experiments. At the same time 
the understanding of the flux and spectrum of Pg from 
power reactors has reached substantial sophistication. 
The first two "long baseline" experiments, Chooz and 
Palo Verde, have amply demonstrated the capability 
of this new breed of detectors, while providing solid ev- 
idence that i>c — v^ is not the dominant channel in the 
atmospheric neutrino oscillation. Kamland, scheduled 
to begin data-taking shortly, will extend the reach for 
small mass-differences to unprecedented levels. Its size 
and background will move reactor-based experiments to 
a new dimension, with several new physics opportunities 
in the essentially background-less detection of v^ from a 
number of natural sources. 
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